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Theoretkal  Studies  of  the  Electronic  and  Resonance 
Structure  of  Atomic  and  Molecular  Negative  Ions 


ABSTRACT 

This  technical  program  constitutes  a  theoretical  research  investigation  of  the  electronic 
structure  of  atomic  and  molecular  ions  and  of  radiative  and  collisional  processes  describing  their 
formation  and/or  destruction.  This  study  is  directed  toward  elucidating  the  structure  of  both 
thermodynamically  stable  cations  and  anions,  and  resonance  or  metastable  excited  states  in  such 
systems.  In  addition,  this  study  will  include  an  analysis  of  the  effects  of  both  static  E-fields  and 
time  dependent  laser  fields  on  resonance  excitation,  stimulated  autodetachment  and  radiative 
electron  attachment  processes  in  negative  ion  systems  and  ion-pair  production  in  Rydberg  atom 
collisions.  The  goals  of  this  research  program  are  to  develop  accurate  computational  techniques 
which  are  applicable  to  studies  of  the  electronic  structure  and  radiative  properties  of  small  atomic 
and  molecular  ions  and  to  carry  out  detailed  calculations  of  the  structure  of  several  prototype 
species,  including  light  element  alkalis,  alkaline  earths  and  the  noble  gases.  It  is  anticipated  that 
the  results  of  this  program  wiU  furnish  fundamental  data  and  provide  a  better  understanding  of 
several  detailed  experimental  studies  that  are  currently  in  progress  in  this  area,  including  work  at 
SRI  and  ORNL.  Areas  of  application  of  this  research  work  include  thermonuclear  reactor 
programs,  neutral  and  negative  ion  particle  beam  source  development,  high  current  switching 
devices  and  discharge  plasma  devices. 
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STATEMENT  Of  WORK 
for 

Contract  F49620-81>-C-0019 

The  Contractor  shaJl  conduct  a  theoretical  research  investigation  of  the  electronic  structure 
of  atomic  and  molecular  ions  and  of  the  radiative  and  collisional  mechanisms  describing  their 
formation  and/or  destruction.  This  investigation  shall  include  ab  initio  calculations  of  atomic  and 
molecular  electronic  structures,  an  analysis  of  resonance  excitation  and  metastable  states, 
stimulated  autodetachment  and  radiative  electron  attachment  processes  in  prototype  ion  systems, 
and  studies  of  ion-pair  production  from  the  collisions  of  Rydberg  atoms.  Experimental  data 
collected  at  ORNL,  LAI^  SRI  and  other  sites  where  there  are  on-going  studies  of  negative  ion 
sources,  shall  be  iteratively  examined  and  integrated  with  the  results  generated  by  this  theoretical 
investigation.  This  investigation  shall  be  carried  out  with  ab  initio  methods  using  electronic 
structure,  radiative  property  and  collisional  cross-section  codes  that  have  previously  been 
developed  at  this  Center,  including  such  modiGcations  as  are  necessary  to  implement  laser  Geld 
effects.  SpeciGcally,  this  investigation  shall  comprise  the  following: 

1.  Quantum  mechanical  calculations  of  the  electronic  structure,  potential  energy  surfaces 
(for  molecular  species)  and  resonance  positions  and  widths  for  positive  and  negative 
atomic  and  molecular  ions.  These  calculations  shall  be  carried  out  for  species  to  be 
chosen  from  the  Grst  and  second  row  atomic  ions.  Heavier  alkali  ions  may  also  be  in¬ 
cluded.  The  choice  of  species  to  be  studied  shall  be  directed  in  part  by  on-going  paral¬ 
lel  experimental  efforts  at  ORNL,  LANL,  SRI,  as  well  as  several  university 
laboratories. 

2.  Variational  R-matrix  calculations  of  resonance  states  excited  through  electron-mole¬ 
cule  collisions. 

3.  Quantum  mechanical  calculations  of  resonance  excitation  and  analysis  of  laser  stimu¬ 
lated  electron  attachment  and  detachment  processes  in  negative  ion  systems.  These 
calculations  shall  include  studies  of  Li~,  Na~  and  He2  ■ 

4.  Analysis  of  ion-pair  formation  in  thermal  energy  Rydberg  atom  coUisions.  Both  alkaU 
and  alkali/noble  gas  systems  shall  be  examined  to  provide  interpretations  of  the  exper¬ 
imental  studies  at  SRI  and  AVCQ  on  these  systems. 
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BACKGROUND 

Recent  successes  in  the  tokamak  program  and  in  other  controlled  thermonuclear  research 
programs  (Ref.  1)  have  focused  attention  on  the  problem  of  developing  an  efficient  high-energy 
particle  beam  source.  For  several  applications,  a  neutral  beam  with  energies  above  200  meV  is 
desired.  The  acceleration  of  negative  ions  to  such  energies,  followed  by  neutralization  through  a 
stripping  reaction,  appears  at  the  present  time  to  be  the  most  efficient  approach  for  producing  a 
high  energy  neutral  particle  beam. 

In  another  area,  the  feasibility  of  particle  beams  (both  charged  and  neutral)  as  military 
defense  weapons  has  been  under  study.  The  proposed  endoatmospheric  applications  require  high 
intensity  sources  and  mainly  focus  on  the  problems  of  beam  stability  and  propagation  characteris¬ 
tics.  Exoatmospheric  applications  require  lower  intensity  sources  but  of  very  high  quality.  Design 
goals  are  highly  collimated  beams  with  a  narrow  energy  spread. 

A  magnetron-type  negative  H"  source  has  been  reported  by  Belchenko,  Dimov  and 
Dudnikov  (BDD)  (Ref.  2)  that  has  produced  H”  current  densities  of  several  A  cm"^  .  This  device 
operates  as  a  plasma  discharge  in  an  atmosphere  of  hydrogen  gas  with  cesium  or  other  alkalis 
present  at  ~  0.01  percent.  The  mechanism  for  the  production  of  H'  was  believed  to  involve  a 
catalyzed  surface  reaction  whereby  H"  ions  are  produced  by  transfer  of  electrons  from  Cs  atoms 
that  are  adsorbed  on  the  cathode  surface  (Refs.  3,  4).  The  detailed  kinetic  mechanisms  of  such 
surface  reactions  are  still  uncertain  and  parametric  experimental  studies  are  currently  underway 
at  LANL,  LBL,  IRT  and  at  other  laboratories  to  elucidate  the  mechanisms  and  operating 
characteristics  of  BDD  and  similar  devices  and  to  develop  information  for  their  scale-up  to 
higher  current  densities.  Diagnostics  of  H”  source  devices  are  also  underway  at  Brookhaven 
(Ref.  5)  using  beam  probe  and  spectroscopic  techniques. 

Concurrent  with  these  surface-plasma  reactions  are  several  electron-  H2  /H  volume- 
dependent  processes  that  may  lead  either  to  the  production  of  H'  or,  in  reverse,  may  act  as 
important  destructive  processes  of  the  negative  ions  after  they  are  formed.  Photodetachment 
experiments  by  Bacal  and  Hamilton  (Ref.  6)  in  hydrogen  plasmas  indicate  H~  densities  100  times 
larger  than  that  predicted  from  simple  electron  attachment  mechanisms.  Further,  these  experi¬ 
ments  indicate  a  nonlinear  dependence  of  the  production  of  H~  on  electron  density,  at  least  for 
densities  less  than  10^*^  cm~^ .  Several  mechanisms  (Ref.  7)  have  been  proposed  to  explain  these 
volume-dependent  H~  production  processes  but  the  cross-sections  for  many  of  the  reactions 
leading  to  H~  are  uncertain. 

In  addition  to  studies  of  thermodynamically  stable  negative  ions,  there  has  been  recent 
interest  in  the  metastable  or  resonant  states  in  these  systems.  Compound  states  of  atomic  negative 
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ions  have  been  a  subject  of  interest  in  atomic  physics  for  over  forty  years.  The  first  observation  of 
a  compound  negative  ion  state  was  the  detection  of  long-lived  He"  ions  in  a  mass  spectrometer 
by  Hiby  (Ref.  8)  in  1939.  Holoien  and  Midtdal  (Ref.  9)  theoretically  showed  that  He"  (Is2s2p/P‘’ 
was  stable  relative  to  He(ls2s)^S  and  thereby  offered  an  explanation  of  Hiby’s  observation  of 
long-lived  He"  ions.  In  1963,  Schulz  (Ref.  10)  reported  evidence  for  a  scattering  resonance  in  He 
at  ~  19.3  eV  which  was  attributed  toHe"(ls2s^)^S .  This  was  the  first  experimental  observation  of 

an  electron  scattering  resonance  in  an  atomic  system.  Since  that  time,  compound  states  of  atoms 
have  been  seen  in  many  elastic  and  inelastic  channels  in  electron-atom  scattering  (Ref.  11),  in 
ion-atom  scattering  cross  sections  (Ref.  12),  in  laser  photodetachment  cross  sections  (Ref.  13), 
and  in  dissociative  electron  attachment  studies  (Refs.  14, 15).  Compound  states  of  doubly  charged 
negative  ions  have  also  been  detected  (Ref.  16). 

The  spectroscopy  of  metastable  negative  ions  has  been  the  subject  of  several  experimental 
and  theoretical  investigations  (Refs.  17-21).  Although  much  attention  has  been  given  to  the 
problem  of  atomic  anion  formation,  more  recent  studies  have  examined  diatomic  anions, 
especially  those  of  the  lighter  elements  (Refs.  22-25).  In  a  recent  study  (Ref.  23),  the  electronic 
structure  of  the  newly  discovered  (Ref.  22)  metastable  He2  ion  was  examined  using  Cl  wavefunc- 

tions.  This  study  indicated  that  only  the  ^  Ilg  state  of  Hez  was  metastable  relative  to 

He2  ]  +  e  .  More  recently,  by  double  charge  exchange  of  He^  ions  in  lithium  vapor,  Kvale, 

et  al.  (Ref.  26)  have  examined  the  autodetachment  spectrum  for  *  He2  and  ^  He2  ions  and  have 
observed  two  autodetachment  peaks.  The  first  (low-energy)  peak  was  identified  as  autodetach¬ 
ment  from  the  first  vibrational  level  of  He2  [^Ilg]  to  He2  [a^2^  ]  +  e  and  the  second  as  autode¬ 
tachment  to  the  He2  +  e  repulsive  continuum,  similar  to  that  observed  by  Bae,  et  al. 

(Ref.  22). 

The  coupling  of  this  experimental  data  to  theoretical  studies  of  the  electronic  structure  of  the 
He2  anion  has  yielded  a  good  understanding  of  the  spectroscopy  of  this  system.  Because  of  the 
difficulties  in  conducting  experiments  to  measure  kinetic  processes  in  these  metastable  and 
compound  ion  systems,  the  development  of  sound  theoretical  methods  for  calculating  the 
positions  and  widths  from  analysis  of  their  electronic  structure  appears  desirable.  In  this 
direction,  we  have  undertaken  a  series  of  calculations  for  the  He/He“/He2  and  Li/Li"/Li2  sys¬ 
tems  in  an  effort  to  understand  their  resonance  structure. 

The  calculations  for  the  He"/He2  system  were  carried  out  to  gain  a  better  understanding  of 
the  ocperiments  of  Bae,  Coggiola  and  Peterson  (Ref.  22),  who  were  the  first  to  report  observations 
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of  the  metastable  He2  ion,  and  the  more  recent  autodetachment  experiments  reported  by  Kvale, 
et  al.  (Ref.  26)  for  this  system.  Our  results  show  that  only  the  lowest  ^  ITg  of  Hea  is  metastable 
relative  to  Hea  [a^Zu  ]  +  e  and  that  the  lowest  lying  ^  !!„,  ^  and  ^  Zj"  states  are  all  autoioniz- 
ing.  However  since  the  spin-orbit  state  has  a  long  (~  400  n  sec)  lifetime,  it  should  be 

possible  to  laser  induce  transitions  from  this  level  to  the  upper  resonant  negative  ion  states.  Such 
an  experiment,  using  a  tunable  laser,  is  now  under  development  by  Compton  and  co-workers  at 
Oak  Ridge  National  Laboratory.  A  complimentary  experiment  would  involve  the  possible 
production  of  thermodynamically  stable  negative  ions  by  stimulated  radiative  attachment  of 
electrons  to  neutral  atoms  or  molecules.  An  alkali  system  would  appear  to  be  a  good  choice  for 
such  studies  since  both  Li~  and  Lia  are  stable  and  exhibit  an  electronic  structure  which  includes 
resonances  that  could  be  reached  by  conventional  tunable  lasers. 

In  a  related  area,  there  is  now  considerable  interest  in  understanding  ion-pair  production  in 
low  energy  collisions  of  Rydberg  atoms.  Ciocca,  et  al.  (Ref.  27)  and  Cheret  and  Barbier  (Ref.  28) 
have  reported  efficient  ion-pair  formation  in  laser  excited  alkali  systems.  In  contrast,  McGeoch 
and  Schlier  (Ref.  29)  find  no  ion-pair  production,  as  a  primary  process,  in  thermal  energy 
collisions  between  lithium  Rydberg  atoms.  Conventional  curve-crossing  mechanisms  are  not 
capable  of  explaining  these  results  and  new  theoretical  developments  are  needed.  An  experiment 
to  study  ion-pair  formation  in  the  reaction  of  He(^S)  atoms  and  alkali  Itydberg  atoms  has  been 

suggested  by  Peterson  (SRI).  This  system  would  be  an  ideal  candidate  for  a  collaborative 
experimental  and  theoretical  study  of  the  mechanism  of  ion-pair  production. 

The  goals  of  this  research  program  were  to  identify  the  electronic  and  resonance  structure  of 
atomic  and  molecular  positive  and  negative  ions,  to  study  the  effects  of  high  laser  fields  on 
stimulated  autodetachment  and  radiative  electron  attachment  in  such  systems,  and  to  carry  out 
theoretical  studies  of  atomic  Rydberg  state  interactions  leading  to  ion-pair  production.  This 
theoretical  research  effort  was  developed  to  compliment  the  experimental  studies  of  ion  formation 
processes  that  are  being  carried  out  by  Peterson  and  co-workers  (SRI),  by  McGeoch  (AVCO)  and 
by  several  European  research  groups. 
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RESEARCH  PROGRESS 

1 .  Electronic  Structure  and  Autodetachment  Studies  of  Molecular  Negative  Ions 

During  the  past  several  years,  we  have  undertaken  a  series  of  calculations  of  Li2  in  order  to 
assess  the  energetics  of  this  anion.  McGeoch  (Ref.  30)  has  suggested  that  a  Li"  source  should  be 
possible  through  a  series  of  reactions  similar  to  those  found  for  H2  .  The  initial  step  involves  the 
formation  of  the  Li2  dimer  and  excitation  to  low-lying  excited  states  via  fast  electron  collisions  or 
photon  pumping.  The  latter  process  is  that  currently  being  investigated  by  McGeoch.  The 
reactions  are: 

e  (fast)  +  Li2  [X  -  U*  [A  *Z+,B  ‘nj  +  e  (1) 

or 

Li2  [X  *Zg+]  +  hv  -  U*  [A  B  ^Uu]  (2) 

The  A  'Zu  and  B  ^Ilu  states  both  exhibit  fast  radiative  decay  back  to  the  ground  X  ^Z^  state 
but  the  Franck-Condon  transitions  favor  a  shift  in  the  population  of  Li2  to  higher  vibrational 
states: 

Li2[A  ^Z+,B  -  U|  pc  ^Zg+  (v  >  5)]  +  hv  (3) 

Our  present  understanding  of  the  energetics  of  the  Li2/Li2  system  is  shown  in  Figs.  1  and  2. 
The  dashed  curves  in  Fig.  1  were  estimates  based  on  H2  and  the  known  molecular  orbital 
structure  of  Li2  •  Our  uiost  recent  ab  initio  calculations,  which  indicate  that  the  lowest  state 

of  Li2 .  is  slightly  attractive  all  the  way  into  the  curve-crossing  with  the  ground  state  of  Li2 .  are 
shown  in  Fig.  2.  This  suggests  that  electron  dissociative-attachment  is  possible  for  this  system  for 
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vibrationally  excited  U2  molecules; 


e  (slow)  +  Li^  [X  —  Li  [^S]  +  Li-f^S] 


(4) 


The  final  potential  energy  curves  governing  reaction  (4)  are  shown  in  Fig.  2.  There  is  still 
considerable  uncertainty  in  the  location  of  the  excited  reasonant  states  of  Li2 ,  especially  in  the 

auto-ionizing  region  above  Li2  +  e.  The  ground  [X  ]  state  is  well  characterized  theoretically 
(Refs.  31,  32)  but  little  is  known  about  the  location  of  the  low-lying  excited  electronic  states  which 
connect  to  Li  [^P]  +  Li“[^S] .  In  addition,  by  analogy  with  H2 ,  low-lying  electronic  states  of  Lii 

with  MO  configuration  (lo|  loj  2ogl;r2)[^’^S^  ]  may,  in  fact,  actually  lie  lower  than  some  of  the 

states  shown  in  Fig.  1.  Such  states  dissociate  to  very  high  atom-ion  limits  and  their  general 
characteristics  are  uncertain.  The  curves  shown  in  Fig.  2  represent  our  best  estimates  at  the 
present  time  but  more  definitive  calculations  are  required  to  define  the  curve-crossing  regions  for 
dissociative-attachment. 

The  lowest  state  of  Li2  has  the  molecular  orbital  configuration  [lo^  2o^  2(7g]  for 
large  separations.  At  short  intemuclear  separations,  this  state  mixes  with  a  second  bound 
state  arising  from  Li  [^P]  +  Li"[*S]  with  a  molecular  orbital  configuration  [lOg  lo^  Tog  3og] . 
This  mixing  at  intermediate  separations  results  in  a  lowering  of  the  interaction  potential  for  the 
lowest  ^Zg  state.  Our  preliminary  calculations,  as  shown  in  Fig.  1,  indicate  that  this  mixing  is 

stronger  than  expected  from  MO  considerations  and  lowers  the  ^Zg  potential  curve  into  a 
feasible  region  for  dissociative-attachment  to  occur. 

The  resonance  structure  associated  with  electron  scattering  by  the  hydrogen  molecule  has 
been  the  subject  of  many  theoretical  and  experimental  studies  (Refs.  11,  33,  34).  In  contrast, 
electron  scattering  by  IJ2  has  received  relatively  little  attention.  Miller,  et  al.  (Ref.  35)  have 
measured  total  cross  sections  for  e  +  Li2  scattering  in  the  0.5-10.0  eV  region  and  find  a  very 

large  (~  500  )  and  rising  total  cross  section  for  the  low-energy  region.  These  results  suggest 

that  some  sort  of  shape  resonance  may  be  found  in  this  system  corresponding  to  a  low-lying  state 
of  Li2 .  Such  resonant  states  have  been  observed  in  electron  impact  excitation  of  Na2  (Ref.  36). 
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In  order  to  clarify  the  electronic  structure  of  the  Li2  anion,  a  series  of  detailed  ab  initio 
calculations  of  the  parent  state  of  Li2  and  of  several  possible  symmetries  of  the  anion  were 
undertaken.  The  calculational  framework  that  was  chosen  was  a  configuration-interaction  (Cl) 
expansion  built  from  orthonormal  Slater-type  orbitals  (STO's).  An  extended  (14cr,  6;r,  2<5)  basis 
was  optimized  using  several  methods  of  approach.  Complete  orbital  optimization  was  possible 
only  for  Li2  and  the  ground  X  ^2,?^  state  of  Li2  since  the  excited  anion  states  are  either 
autodetaching  at  some  critical  intemuclear  separation  or  are  pure  resonant  states.  The  lowest 
^2^  state  of  Li2 ,  which  is  the  principle  channel  for  DA  of  e  +  Li2 ,  is  variationally  stable  for 
large  intemuclear  separations,  and  a  smooth  extrapolation  of  optimized  orbital  exponents  for  this 
state  was  carried  out  into  the  autodetaching  region.  Hie  purely  resonant  states  of  Li~  and  Li2 
were  treated  using  a  modified  nuclear  charge  hamiltonian,  similar  to  that  described  by  Nestmann, 
et  al.  (Ref.  37).  In  this  method  the  resonant  states  are  variationally  stabilized  by  an  increased 
nuclear  charge  and  an  analytic  extrapolation  is  carried  out  to  determine  the  real  part  of  the 
expectation  value  of  the  unperturbed  hamiltonian. 

Our  calculated  electron  affinity  for  Li2  >  derived  from  our  location  of  the  X  ^2^*  state  of 
Ul ,  is  0.42  ey  in  good  agreement  with  previous  theoretical  studies.  The  A  ^2g  state  of  Liz 

exhibits  a  complex  potential  energy  for  intemuclear  separations  smaller  than  3.45  A ,  where  a 
crossing  with  X  ^2^  state  of  Liz  found.  This  resonant  state  of  Liz  is  mainly  of  the  Feshbach 
type  and  has  a  long  lifetime  against  autoionization.  We  find  a  small  width  for  decay  to  the  ground 
state  of  Liz ,  r  (R)  0.01  k(R).  Of  special  interest  is  the  low-lying  ^Ilu  resonant  state  which  lies 

appr(»dmately  150  meV  above  the  ground  state  of  Liz  •  ^is  anion  state  has  a  dominant  MO 
configuration  of  (lOg  ^  lou  ^2ag  ^  Ltu)  ,  indicating  a  relatively  large  capture  width  for  low-energy 
electron  scattering,  and  is  probably  responsible  for  the  large  scattering  cross  sections  observed  by 
Miller,  et  al.  (Ref.  35)  for  low  collisional  energies.  Hie  character  of  this  ^IIu  state  is  similar  to  the 

^Ilg  resonant  state  of  Nz  in  that  the  anion  state  lies  lower  than  the  corresponding  neutral 
parent  for  a  wide  region  of  intemuclear  separations  (3.5  <  R  <  9.0  A ). 

Further  theoretical  studies  of  these  interaction  potentials  are  required  to  define  the  overall 
reactive  kinetics  leading  to  Li~  production.  A  study  parallel  to  that  reported  by  Hiskes  (Ref.  38) 
for  Hz  has  been  carried  out  to  estimate  the  overall  cross  sections  for  production  of  vibrationally 

excited  Liz  ■  In  this  study,  cross-sections  for  vibrational  recitation  (VE)  of  the  X  ^2^  state  of 
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Li2  via  formation  of  the  excited  A  or  B  *!!„  electronic  states  were  reported.  For  vibration¬ 
al  excitation  through  the  A  *Z,J'  state,  the  cross  sections  are  nearly  constant  for  forming 
3  <  v"  <  9  via  electron  collisional  excitation.  For  photon  pumping  (670  nm  <  A  <  700  nm)  of 

the  A  state,  levels  v"  <  9  are  predominantly  formed.  Vibrational  excitation  via  the  B  ^Ilu 
state  has  a  lower  probability.  In  addition,  theoretical  studies  are  in  progress  to  estimate  the 
low-energy  electron  dissociative  attachment  cross-section  for  e  +  Li^  Li  +  Li~  .  Our  initial 
results  are  shown  in  Fig.  3  which  indicate  an  attachment  cross  section  >  10"^^  cm^  for  vibration- 
ally  excited  Li2  •  The  enhancement  in  the  attachment  cross  section  with  vibrational  energy  is 
similar  to  that  found  for  the  e  -h  H2  attachment. 

In  addition  to  our  studies  of  the  Li2/Li2  system,  additional  calculations  were  carried  out  on 
the  He2  metastable  anion.  The  existence  of  metastable  negative  ions  has  been  the  subject  of 
several  recent  experimental  and  theoretical  investigations  (Refs.  19-21).  Although  much  attention 
has  been  given  to  the  problem  of  atomic  anion  formation,  more  recent  studies  have  examined 
diatomic  anions,  especially  those  of  the  lighter  elements  (Refs.  22-25).  In  a  recent  study  (Ref.  23), 
we  have  examined  the  electronic  structure  of  the  newly  discovered  (Ref.  22)  metastable  He^  ion 

using  Cl  wavefunctions.  This  study  indicated  that  only  the  ^Flg  state  of  the  He2  was  metastable 
relative  to  He2  [a^Zj"  ]  +  e ,  with  an  indicated  electron  affinity  of  0.233  eV  and  a  vibrational 
frequency  (cUe)  of  1783.0  cm”^ .  A  low-lying  ^Zg  resonance  state  of  Hei  was  also  found  but  it 

was  not  possible,  within  the  error  of  the  calculations,  to  ascertain  whether  or  not  this  state  might 
also  be  metastable. 

More  recently,  by  double  charge  exchange  of  He^  ions  in  lithium  vapor,  Kvale,  et  al.  (Ref. 

26)  have  examined  the  autodetachment  spectrum  for  ^He^  and  ^Hei  ions  and  have  observed 
two  autodetachment  peaks.  The  first  (low-energy)  peak  was  identified  as  autodetachment  from 
the  first  vibrational  level  of  He2  [^Ilg]  to  He2  (a^Z,r  ]  +  e  and  the  second  as  autodetachment  to 

the  He2  [^Zg  ]  -I-  e  repulsive  continuum,  similar  to  that  observed  by  Bae,  et  al.  (Ref.  22).  This 
experiment  by  Kvale,  et  al.  suggested  that  either  the  low-energy  autodetachment  was  from  the  first 
excitcu  vibrational  level  of  the  ^Ilg  state  as  postulated,  or  from  a  second  metastable  electronic 

state  (possibly  ^Z^  )  lying  in  this  same  energy  range.  In  either  case,  our  earlier  results,  which 
indicated  that  the  first  excited  vibrational  level  of  the  ^Ilg  state  was  non-autodetaching  and  that 
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the  lowest  *  state  was  purely  resonant  in  character,  were  in  disagreement  with  these  exper¬ 
imental  observations. 

In  order  to  clarify  the  spectroscopy  of  Hea ,  a  re-examination  of  both  the  a  state  of 
Hea  and  the  low-lying  symmetries  of  Hea  was  carried  out.  Our  previous  studies  of  this  system 
were  based  on  Cl  wavefunctions  constructed  from  a  (14<7,  6n,  2&)  STO-MO  basis.  The  atomic 

orbital  parameters  were  taken  from  optimized  calculations  for  the  ^P°  state  of  He” .  In  these  new 

calculations,  a  complete  non-linear  orbital  optimization  was  carried  out  at  the  equilibrium 
intemuclear  separations.  This  non-linear  optimization  was  performed  separately  for  both  the 
neutral  and  negative  ion  states  in  an  attempt  to  quantitatively  characterize  the  small  energy  shifts 
arising  from  polarization  effects  in  the  formation  of  the  anion.  This  approach  permitted  examina¬ 
tion  of  smaller  changes  in  the  charge  distribution,  upon  anion  formation,  than  would  have  been 
possible  by  fixed  basis  set  expansions  with  exponents  chosen  to  insure  freedom  from  linear 
dependency  problems. 

For  studies  of  the  lowest  ^Z^  state,  a  modified  nuclear  charge  hamiltonian  was  employed 
similar  to  the  method  described  by  Nestmann  and  Peyerimhoff  (Ref.  39).  Stabilization  of  this  state 
(relative  to  He^  (a^Z,f  ]  +  e )  was  realized  by  increasing  the  nuclear  charge  by  ~  0.1%.  A  careful 
study  of  the  non-linear  parameters  as  a  function  of  increased  nuclear  charge  clearly  indicated  that 
the  lowest  ^Z^  state  is  purely  resonant  in  character  and  that  this  state  is  autoionizing  relative  to 

He2  [a^Z^  ]  +  e .  We  find  that  this  lowest  ^Z^  state  lies  0.074  eV  above  He2  [a^Zu  ]  and  has  a 
dominant  MO  structure  of  (lso|  Isou  2sag  2sau) .  A  preliminary  analysis  of  the  complex  part  of 

the  potential  indicates  that  this  state  has  a  large  width  for  decay  to  the  He2  [a^ZiT]  -i-  e 
continuum  and  is  unlikely  to  be  observable  on  the  time  scale  of  the  detection  apparatus  in  the 
double  charge  exchange  experiments.  Further  calculations  on  the  ^Zg  and  other  low  lying 
resonances  in  He2  are  required  to  define  the  real  and  complex  behavior  of  these  states. 

The  lowest  ^Ilg  state  of  He2  was  found  to  be  variationally  stable  for  all  intemuclear 
separations  that  were  examined.  Our  final  potential  curve  for  this  state  is  displayed  in  Fig.  4.  We 
find  an  electron  affinity  of  0.182  eV  for  the  ^Ilg  state  of  He2 ,  a  value  ~  50  meV  smaller  than  in 

our  previous  calculations  which  employed  an  STO  basis  that  was  optimized  for  the  ^  P°  (Is2s2p) 

state  of  He” .  We  attribute  this  smaller  calculated  electron  affinity  to  our  observation  that  the  use 
of  optimized  basis  functions  resulted  in  a  larger  improvement  in  the  calculated  energy  of 
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He2  {a^Su  ]  than  to  the  corresponding  Hei  [^Ilg]  anion.  Our  predicted  electron  affinity,  coupled 

with  a  calculated  vibrational  spacing  of  1702  cm'^ ,  now  indicates  that  the  first  excited  vibrational 
level  of  the  ^Ilg  state  of  Hea  is  autoionizing  and  is  probably  responsible  for  the  observed  low 
energy  peak  in  the  electron  detachment  studies  of  Kvale,  et  al.  (Ref.  26).  Their  experimental 
estimate  (0.175  eV)  of  the  electron  affinity  of  the  ^Ilg  state  of  He2 ,  relative  to  He2  ] ,  and 

our  calculated  value  (0.182  eV)  are  within  the  uncertainties  of  both  the  experimental  measure¬ 
ments  and  calculations. 

In  a  related  study  at  SRI,  Bae,  Coggiola  and  Peterson  (Ref.  22)  have  measured  the  energy 
spectrum  resulting  from  electron  autodetachment  from  the  lowest  ^Ilg  state  of  He2 .  This 
continuum  sfiectrum  exhibits  a  single  peak  at  15.7  eV  resulting  from  a  vertical  transition  from 
He2  to  the  repulsive  ground  state  of  He2 .  Both  the  shape  and  position  of  this  autodetachment 
spectrum  were  verified  in  later  studies  reported  by  Kvale,  et  al.  (Ref.  26).  However,  neither  the 
shape  or  position  of  the  experimental  spectrum  seemed  to  fit  predictions  based  on  our  He2 
calculations  and  the  generally  accepted  ground  state  He2  potential  reported  by  Foreman,  et  al. 
(Ref.  40). 

In  order  to  resolve  this  discrepancy,  Michels  and  Peterson  (Ref.  41)  jointly  analyzed  the  He2 
continuum  spectra  reported  by  Huffman  et  al.  (Ref.  42)  and  concluded  that  the  Foreman  ground 
state  potential  for  He2  was  essentially  correct.  An  exhaustive  study  by  Michels  and  Peterson  (Ref. 
43)  of  possible  shapes  of  the  potential  energy  curves  for  Hei  resulted  in  the  conclusion  that  the 
He2  beam  must  have  been  rotationally  hot  in  the  SRI  experiments.  A  suitable  fit  to  the  exper¬ 
imental  spectrum  was  obtained  by  assuming  a  rotational  temperature  of  ~  15000  K,  with 
significant  contributions  fi'om  higher  vibrational  levels  of  He2 .  This  collaborative  experimental 
and  theoretical  effort  solved  a  formidable  problem  but  introduced  new  questions  concerning  our 
understanding  of  electron  autodetachment  processes.  Our  conclusions  that  vibrational  autode¬ 
tachment  must  be  slow  in  the  He2  system  runs  counter  to  previous  ideas  and  suggests  that 
electron  coupling  mechanisms  must  be  reexamined. 

2.  Hydrogen  Cation  and  Anion  Cluster  Formation  Studies 

There  have  been  numerous  thermochemical  and  kinetics  studies  of  weakly  bound  ionic 
clusters  during  the  past  decade  (Refs.  44  and  45).  Of  particular  interest  are  the  cation  and 
Hq  anion  sequences  since  they  represent  examples  of  the  chemistry  of  strong  chemical  bonding, 
in  ,  the  simplest  hydrogen  bonded  molecular  ion,  in  ,  and  the  simplest  atomic  anion. 
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H'  .  In  addition,  for  the  first  few  members  of  these  sequences,  detailed  ab  initio  quantum 
mechanical  calculations  of  their  equilibrium  bond  lengths  and  vibrational  spectra  can  be  carried 
out. 


Very  recently,  the  possibility  of  forming  large  HjJ"  clusters  (n  >  1000),  that  could  be  kept 
for  long  periods  of  time  in  a  Penning  trap,  has  been  discussed  (Ref.  46).  Such  a  system  is  a 
prototype  for  bulk  antiproton  or  antihydrogen  storage,  with  uses  as  high  energy  density  material 
or  as  an  energy  source  for  advanced  propulsion  concepts.  The  mechanism  for  formation  and 
storage  of  these  large  cluster  ions  can  be  understood  only  if  the  thermodynamics  and  kinetics  of 
the  cluster  growth  processes  are  identified  and  analyzed  in  detail.  The  purpose  of  our  study  was 
to  examine  the  thermodynamic  stability  of  the  first  few  members  of  the  hydrogen  cation  and 
anion  cluster  sequences  since  their  stability  is  a  requirement  for  defining  a  useful  association 
pathway  for  the  growth  of  large  clusters. 

Numerous  experiments  (Refs.  47  and  48)  have  measured  the  stability  of  the  odd-membered 
hydrogen  clusters  that  are  formed  by  the  association  reaction: 


H+  +  H2  -  H„+^.2  ;  n  =  3,  5,  7,  .  .  .  (5) 

These  clusters  exhibit  stabilization  energies  of  approximately  3-6  kcal/mol,  in  good  agreement 
with  theoretical  calculations  (Refs.  49  and  50).  Until  recently,  the  even-membered  hydrogen 
clusters  have  not  been  observed,  with  the  exception  of  the  suggestion  of  a  stable  H4  ion  formed 
by  the  attachment  of  a  H  atom  to  the  H3  ion  (Ref.  51).  A  very  recent  study  of  the  even-mem¬ 
bered  clusters  (Ref.  52)  suggests  that  the  observation  could  be  interpreted  as  H2D''’ .  This 
study  also  indicates  a  high  degree  of  stability  for  . 

In  contrast,  the  negative  hydrogen  ion  clusters  have  never  been  observed  and  a  recent  careful 
search  for  stable  or  metastable  H2  and  H3  produced  negative  results  (Ref.  53).  There  have  been 
several  theoretical  studies  (Refs.  54  and  55)  which  support  the  existence  of  stable  H3 ,  and  higher 
clusters,  but  more  accurate  <16  initio  calculations  (Ref.  56)  shed  doubt  on  this  conclusion. 

In  order  to  establish  the  thermodynamic  stability  of  these  hydrogen  cluster  ions,  ab  initio 
calculations  were  carried  out  to  detenm'ne  their  electronic  structure  and  stable  geometries.  Several 
Gaussian  basis  sets,  with  increasing  flexibility,  were  employed  in  these  studies  and  levels  of  theoiy 
up  to  single  and  double  configuration  interaction  expansions  (CISD)  were  examined. 
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The  results  of  these  calculations  are  summarized  in  Table  1,  which  includes  the  H”  and  H2 
constituents  as  well  as  H3  .  It  is  clear  from  Table  1  that  the  inclusion  of  the  vibrational  zero  point 
energy  is  necessary  to  assess  the  stability  of  H3 .  We  find  with  every  basis  and  level  of  theory  that 
the  H3  ion  is  not  stable  with  respect  to  H’  and  H2 .  This  conclusion  is  supported  by  the  results 
of  earlier  theoretical  work  (Refs.  55  and  56),  and  by  the  recent  experimental  study  of  Bae, 
Coggiola,  and  Peterson  (Ref.  53). 

Our  calculations  appear  to  rule  out  the  formation  of  hydrogen  anion  clusters  through  either 
of  the  association  reactions: 


Hi  +  H  -  (6) 

+  H2  -  H;;+2  (7) 


Although  previous  studies  by  Hirao  and  Yamabe  (Ref.  56)  indicate  a  weak  stability  for  Hn 
clusters,  up  to  n  =  13,  the  non-existence  of  H2  and  H3  precludes  formation  of  higher  clusters 
through  eit  er  reaction  (6)  or  (7).  There  remains  the  possibility  (albeit  remote)  that  association  of 
H~  with  a  stable  Rydberg  state  of  H3  could  produce  H4 . 

3.  Resonant  Electron  Structure  of  Molecular  Anions 

Electronic  structure  calculations  of  resonant  or  autoionizing  states  are  usually  carried  out  by 
stabilization  methods  or  by  the  use  of  modified  hamiltonians  which  introduce  complex  coordi¬ 
nates.  The  latter  methods  have  the  advantage  that,  for  molecular  systems,  both  the  real  and 
complex  parts  of  the  potential  energy  surface  can  be  directly  defined.  The  disadvantages  include 
difficulty  of  application  and  for  some  systems,  such  as  Li~  and  N2 ,  an  ambiguity  in  the  results  of 
several  approaches.  In  a  recent  paper  (Ref.  39X  Nestmaim  and  Peyerimhoff  have  developed  a 
modified  hamiltonian  method  that  can  be  applied  to  low-lying  resonances,  basically  within  a 
standard  Cl  computational  fi-amework.  In  this  method,  stabilization  of  a  negative  ion  state  is 
achieved  through  a  modification  of  the  nuclear  charge,  z  =  z(l  +  A) .  For  a  critical  value  of 

A  -  Ac ,  a  particular  resonance  can  be  stabilized  below  the  autoionization  threshold.  The 
quantity 


E(A)  =  <  V's  W-l)!  %  > 


(8) 
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can  be  shown  to  be  a  quadratic  function  of  A  and  a  series  of  calculations  for  A  >  Ac  yields 
successive  approximations  to  E  (A)  that  can  be  extrapolated  to  yield  Eq  =  E  (A  =  0) ,  which  is 
the  real  part  of  the  resonance  energy.  Development  of  the  complex  component  can  be  accom¬ 
plished  by  a  partitioning  technique  whereby  the  total  space  of  the  system,  {V'i}  i  =  1.  n,  is 

decomposed  into  a  P  space  spanned  by  the  functions  which  describe  the  neutral  target  and  a  free 
electron  and  a  Q  space  that  constitutes  the  compact  component  of  the  resonance.  The  coupling, 
kj ,  between  the  compact  component  of  the  resonance  and  the  scattering  states  described  with  the 
subspace  P  can  be  written  as: 


kftv  “ 


n-l  1/2  n-1 

Xk,V(n-l)  -  2;‘=f(Ei-Eo)V(n-l) 

j=i  J  |_i=i 


1/2 


(9) 


An  appropriate  resonance  width  can  then  be  obtained  from  the  Fermi  golden  rule 


r(e)  =  fer^kiv  I2AE]1/2 


(10) 


where  AE  is  the  difference  between  the  resonance  state  energy  and  that  of  the  neutral  target. 

We  have  applied  this  formalism  to  the  low-lying  resonances  in  the  He/He'/He^  and 
Li/Li2/Li2  systems  with  good  success.  Additional  studies  of  the  He2  resonances  are  needed 
since  this  system  is  a  candidate  for  studying  laser  induced  transitions  from  a  metastable  level. 
Other  ions  include  Be"  where  the  (ls^2s2p^  -•  (ls^2s^2p)  transition  is  of  interest,  metastable 
levels  in  the  alkaline  earths  and  molecular  anions  such  as  HeNe",  HeH"  and  high-lying  reso¬ 
nances  in  H2 .  Several  of  these  systems  are  currently  being  studied  experimentally  by  Compton 
and  co-workers  at  ORNL,  and  by  Peterson  and  co-workers  at  SRI,  by  means  of  autodetachment 
experiments. 

4.  R-matrix  Calculations  of  Elactron-Molscula  Electronic  Excitation  Collisions 

Our  studies  of  vibrational  excitation  of  IJ2  •  via  electron  or  photon  pumping  of  the  electroni 

cally  excited  states,  A  ^  and  B  ,  resulted  in  relative  cross  sections,  since  the  electronic 

excitation  cross  sections  are  unknown  for  this  system.  Although  excitation  of  the  electronic  states 
of  a  molecule  by  electron  impact  is  one  of  the  simplest  basic  processes  in  molecular  collision 
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phenomena,  progress  in  the  theory  of  this  problem  has  been  slow.  For  a  few  diatomic  molecules, 
systematic  studies  (Refs.  57-64)  of  the  excitation  cross  sections  have  been  made  by  means  of  the 
Born  or  modifled  Born  approximations.  Other  approaches  have  been  the  R-matrix  methods  (Refs. 
65  and  66),  distorted-wave  approximation  (Ref.  67),  the  T-matrix  method  (Ref.  68),  and  the 
close-coupling  methods  as  applied  to  e  +  H2  collisions  by  Chung  and  Lin  (Ref.  69). 

For  electronic  excitation  of  a  molecule  through  collisions  with  electrons,  the  R-matrix 
method  and/or  close-coupling  techniques  would  appear  to  be  the  most  fruitful  approaches  for 
low-energy  collisions  (E  s  5  eV).  The  goal  of  this  research  task  is  to  carry  out  e  +  H2  and  e  + 
Li2  calculations  to  effective  convergence  with  respect  to  both  the  orbital  basis  set  and  the 
representation  of  the  target  molecule  and  scattering  wave  functions.  A  primary  goal  is  the 
calculation  of  electronic  excitation  cross-sections  to  the  low  lying  and  ^Ilu  states.  Con¬ 
verged,  fixed-nuclei,  electron  scattering  results  of  this  kind  are  required  for  subsequent  calcula¬ 
tions  of  nonadiabatic  effects  in  vibrational  excitation  and  dissociative  attachment. 

Computer  programs  designed  for  efficient  quantitative  calculations  of  electron-molecule 
scattering  cross  sections  have  been  constructed  by  scientific  collaborations  involving  several 
research  groups.  The  Belfast/Daresbuiy  group  (Refs.  66,  70)  has  modified  the  ALCHEMY  system 
of  bound  state  electronic  wavefunction  programs  for  molecular  R-matrix  calculations.  Results 
obtained  with  Slater  orbitals  and  an  artificial  center  within  the  molecule  (Ref.  71)  indicate  that 
this  approach  is  feasible  only  for  low  impact  energies,  because  of  slow  convergence.  Oscillatory 
continuum  basis  orbital  functions  are  needed  to  improve  convergence.  In  calculations  of  e  +  N2 
scattering,  (Ref.  72)  numerical  solutions  of  a  spherical  model  potential  scattering  problem  were 
used  as  basis  orbitals,  with  satisfactory  results  up  to  intermediate  energies. 

In  collaboration  with  R.  K.  Nesbet  of  IBM,  we  have  begun  an  extension  and  modification  of 
the  present  R-matrix  code  to  the  CRAY  class  of  computers.  In  this  project,  the  original  IBM- 
based  code  will  be  converted  to  the  CRAY  and  eventually  expanded  to  include  electron-polyatom¬ 
ic  molecule  scattering  within  the  GAUSSIAN  86  framework. 

The  present  program  uses  two  forms  of  continuum  basis  orbitals.  Numerical  asymptotic 
functions  (NAFs)  are  generated  by  solving  the  coupled  multichannel  potential  problem  defined  by 
the  static  and  transition  multipole  moments  of  a  chosen  set  of  target  molecular  states  and 
pscudo«tates.  V^ational  R-matrix  calculations  include  two  independent  NAFs  for  each  open 
scattering  channel,  as  energy-dependent  basis  functions.  Effective  completeness  within  the 
R-matrix  boundary  ri  is  achieved  by  including  spherical  Bessel  functions  that  vanish  at  ri , 

which  by  themselves  define  a  complete  set  for  each  orbital  /-value.  NAFs  are  integrated  inwards 
to  some  inner  point,  ro ,  inside  the  target  charge  distribution,  where  they  are  matched  to  regular 
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spherical  Bessel  functions.  With  two  NAFs  per  open  channel,  the  exact  multichannel  solution  of 
the  asymptotic  scattering  problem  is  brought  inside  the  R-matrix  boundary.  Hence  the  energy-in¬ 
dependent  continuum  basis  functions  (Bessel  functions)  are  needed  only  to  represent  short-range 
and  nonlocal  effects. 

Programs  are  under  development  to  compute  static  and  transition  multipole  moments  and 
the  corresponding  multipole  pseudostates.  Pseudostate  orbitals  are  used  to  describe  the  dynami¬ 
cal  polarization  response  of  the  target  molecular  states.  New  techniques  due  to  Gailitis  (Ref.  73) 
are  used  to  integrate  the  asymptotic  equations. 

Applications  of  this  new  code  to  inelastic  e  -I-  H2  scattering  was  carried  out  first  to  check 
the  convergence  rate  for  the  wavefunction  expansion.  Comparisons  with  the  theoretical  work  of 
Baluja  et  al.  (Ref.  74)  and  Collins  and  Schneider  (Ref.  75)  was  made  along  with  analysis  of  the 
experimental  excitation  cross  section  data.  Attention  will  then  be  focused  on  the  e  +  Lia  system. 
In  particular,  cross  sections  was  calculated  for  vibrational  and  electronic  excitation.  For  low 
energy  electrons,  direct  vibrational  excitation  probably  occurs  through  formation  of  a  resonant 
Li  2  structure; 


e  +  y  (v")  -*  Ua  -  e  Lif  (v  >  v")  (11) 

High  energy  electrons  excite  the  vibrational  levels  of  Lia  via  the  indirect  mechanism  that  we  have 
previously  described  (Ref.  76): 

e  Ua  [X  ‘Zg+J  -  y*  [A  B  ’Hu]  +  e  (12) 

Ul  [A  B  -*  y|  [X  ^Zg^]  -t-  hv  (13) 

For  a  completely  ab  initio  analysis  of  the  dissociative  attachment  reactions. 


e  +  y|  [x‘Zg^I  -  y  i^s]  +  y-  [*s]  (i4) 

the  excitation  cross-sections  corresponding  to  Eqs.  (IIX  (12)  and  (13)  must  be  calculated.  Further 

applications  include  the  e  +  He*  [a^Z^  ]  system  and,  if  computationally  practical,  the  low-lying 
resonance  structure  of  e  +  HaO  collisions.  This  research  on  dissociative  attachment  reactions  is 
being  carried  out  in  collaboration  with  Prof.  J.  M.  Wadehra  of  Wayne  State  University. 
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5.  Laser  Induced  Transitions  in  Negative  ions 

In  the  presence  of  a  strong  laser  field,  the  collision  of  an  electron  with  a  neutral  atom  (or 
molecule)  that  has  a  thermodynamically  stable  anion  may  result  in  resonant  attachments.  This 
process  was  first  suggested  by  R.  Compton  (Ref.  77)  and  has  been  analyzed  in  some  detail  by 
Mittleman  (Ref.  78).  Using  Li"  as  a  prototype  system  the  overall  process  would  be: 


e  +  Li  [2S(2s)]  -  Li"  [’P(2s2p)] 


(15) 


laser 


E(t)  sin  vt 


+  Li-[*P(2s2p)] 


Li-[*S(2s2)]  +  hv 


(16) 


For  monoenergetic  electrons,  with  Ee  =  E(Li"  *P)  -  E(Li  ^S),  stimulated  emission  from  the 
resonant  *P  state  of  Li"  should  occur  for  a  laser  frequency  such  that 


hv  =  E(U-  ip)  -  E(y-  *S)  =  Ee  +  Ea  . 


(17) 


The  detailed  analysis  of  laser  induced  radiative  attachment  (and  resonance  excitation) 
requires  a  field  dependent  hamiltonian  such  that  relaxation  effects  of  the  neutral  atom  are  treated 
in  a  self-consistent  fashion.  One  approach  that  might  prove  fruitful  is  a  modification  of  a 
variational  scattering  formalism  to  include  induced  field  effects. 

We  have  implemented  an  expansion  technique  for  inelastic  electron-atom  scattering  (Ref.  79) 
that  might  be  applied  to  the  field-dependent  case  provided  certain  bound-free  integrals  can  be 
evaluated.  Our  experience  in  this  area  dates  back  some  twenty  years  and  the  integral  problem  can 
probably  be  solved  in  a  satisfactory  maimer. 

Our  expansion  technique  formalism  follows  closely  along  the  developments  of  Kohn  and 
Hulthen.  For  the  central-field  (laser  off)  Hamiltonian  K  of  finite  range  (as  would  be  encountered 
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in  electron-neutral  particle  scattering)  and  at  a  given  energy  E,  we  consider  approximate  stationa' 
ry-state  wavefunctions  yj ,  of  the  general  form. 


2n  n 

y>  =  Y.  +  X  h\1Ji  (18) 

i=l  i=l 

Here  the  are  the  complete  set  of  asymptotic  eigenfunctions  (unbound  states)  appropriate  to 
energy  E,  and  the  rji  are  quadratically  integrable  short-range  functions  needed  to  approximate  rp 
within  the  range  of  36 .  In  most  implementations  of  the  R-matrix  method,  only  the  set  of  Tfi  are 

included  but  with  the  addition  of  conditioned  functions  capable  of  an  approximate  description  of 
the  asymptotic  region. 

In  a  n-channel  problem  there  will  be  n  different  stationary-state  wavefunctions  to  be 
approximated  by  Eq.  (18).  All  scattering  processes  can  be  characterized  from  the  ai  coefficients 
of  these  wavefunctions  since  the  ai  define  the  relative  amplitudes  and  phases  in  the  various 
channels.  To  variationally  determine  these  wavefunctions,  we  assume  the  expansion 

2n  n 

(3t-E)V'  =  2  ^  2 

i=l  i=l 

Equation  (19)  amounts  to  the  assumption  that  short-range  functions  beyond  those  important  in 
Ip  do  not  contribute  appreciably  to  (36  -E)  y; .  Combining  Eqs.  (18)  and  (19)  and  taking  scalar 

products  with  i/j  and  ^ ,  we  have 


M,^  a  +  M,p,  b  =  d  (20) 

a  +  b  =  Nc 

where  the  Mi^,  +  M ,  etc.,  are  matrices  of  <  7i|36-E|^  >,  <  7i|36-E|j7j  >  ,  etc.  The 
solution  of  Eqs.  (20)  and  (21)  can  be  accompUshed  in  several  approximations  (Refs.  79  and  80). 
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Introduction  of  the  laser  field  as  a  perturbation  to  the  n  + 1  electron  hamiltonian  %  now 
requires  that 


|^i^-M'(hv)j  V'Xt)  =  0  (22) 

with  36'  =  36  -  E(t)  ^  •  Pn  and  subject  to  the  constraint  that  the  central-field  hamiltonian 
satisfy  a  self-consistent  minimization  condition 


(36  -  Ee)  yf  =  minimum 


(23) 


The  function,  ip  =  i  =  1,  .  .  .  n,  can  be  decomposed  into  a  P  space  spanned  by 

functions  which  describe  the  distorted  neutral  target  plus  the  electron  and  a  Q  space  that 
describes  the  compact  resonant  state.  Then,  to  first  order,  the  transition  from  the  resonant  part  of 
KV'res  =  Q  36  Q ,  to  the  final  negative  ion  state  ipt  is  governed  by  the  matrix  element 

<  W |e  •  Pn  IV'res  > 

where  ifH  is  variationally  constructed  from  the  functions  {V'i}  i  =  .  .  .  n  and  projected  to  be 

orthogonal  to  the  resonant  state. 

Although  modified  field  hamiltonians,  within  a  self-consistent  formulation,  have  been  used 
in  calculations  of  dipole  polarizabilities,  the  time  dependance  introduced  through  Eq.  (22)  must 
now  be  treated  by  successive  approximations.  The  matrix  element  calculations  remain  a  formida¬ 
ble  task. 

6.  Ion-Pair  Formation  in  Rydberg  Atom  Coiiiaions 

The  negative  ions  of  light  atoms  are  currently  being  studied  for  their  possible  application  in 
gaseous  discharges,  fusion  plasmas  and  gas  lasers  (Ref.  81).  Our  source  for  the  volume  production 
of  atmnic  anions  is  the  process  of  dissociative  electron  attachment  to  molecules  (Ref.  82).  This 
process  has  been  studied  experimentally  for  low  energy  electron-hydrogen  molecule  collisions  by 
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Allen  and  Wong  (Ref.  83).  A  parallel  theoretical  study  has  been  reported  by  Wadehra  and 
Bardsley  (Ref.  84).  More  recently,  McGeoch  and  Schlier  (Ref.  85)  have  examined  dissociative 
attachment  (DA)  in  electron-lithium  molecule  collisions  and  have  found  large  DA  rates  for 
attachment  to  highly  vibrationally  excited  Li2  molecules.  The  effect  of  vibrational  excitation  on 
the  DA  rates  has  been  studied  by  Hiskes  (Ref.  86)  for  e  +  H2  collisions  and  by  Wadehra  and 
Michels  (Ref.  87)  for  the  e  +  Li2  system. 

In  the  e  +  Li2  system,  it  has  been  shown  (Ref.  87)  that  excitation  of  Li2  to  the  A 
state,  either  by  electron  impact  or  by  photon  pumping,  results  in  an  enhanced  vibrational 
distribution  upon  radiative  decay  of  the  A  state  to  the  ground  X  state  of  Li2  • 

cross  sections  for  vibrational  excitation  of  Li2  X  *2/  via  low  energy  electron  collisional  excita¬ 
tion  of  the  A  ^2,^  state  exhibit  a  relatively  flat  vibrational  distribution  in  the  range  3  ^  v"  ^  9 . 
This  suggests  that  Li2  molecules  which  are  vibrationally  excited  to  the  DA  threshold  for  Li*  for¬ 
mation  should  exhibit  the  largest  attachment  cross  sections,  a  result  that  is  indicated  by  the 
studies  of  McGeoch  and  Schlier  (Ref.  85). 

An  alternative  mechanism  for  negative  ion  formation  was  suggested  by  Lee  and  Mahan  (Ref. 
88),  who  proposed  that  ion-pair  formation: 


X*  X  —  X+  +  X- 


(25) 


should  occur  for  systems  where  the  energetics  for  ion-pair  formation  are  competitive  with 
associative  ionization: 


X*  X  —  X^  +  e. 


(26) 


The  alkali  metals,  including  lithium,  are  therefore  good  candidates  for  such  a  study  of  ion-pair 
formation.  Ciocca,  et  al.  (Ref.  27)  have  recently  reported  on  the  formation  of  Na~  ions  in  highly 
excited  (n  =  7-40)  Rydberg  atom  collisions.  In  a  similar  experiment,  Cheret  and  Barbier  (Ref.  28) 
have  reported  efficient  ion-pair  formation  in  collisions  between  ground  state  and  Rydberg  Rb 
atoms  formed  by  laser  excitation.  A  curve  crossing  mechanism  was  suggested  to  explain  these 
results. 
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More  recently,  McGeoch  and  Schlier  (Ref.  29)  have  looked  for  ion-pair  formation  in 
collisions  between  two  Rydberg  Li  atoms  (n  =  5-9).  The  most  prominent  processed  observed  was 
associative  ionization: 


Li**  +  U**  -  U2^  +  e  (27) 

rather  than  ion-pair  production,  which  was  not  observable  at  a  rate  ~  300  times  less  than 
associative  ionization.  McGeoch  and  Schlier  followed  the  time  dependence  for  the  formation  of 
Li~  and  find  that  it  follows  the  production  route  by  associative  ionization  of  e  +  Li2 .  These 
results,  for  the  first  time,  rule  out  direct  ion-pair  formation  as  an  important  process,  at  least  in  Li 
atom  collisions,  and  suggest  that  the  earlier  work  of  Ciocca,  et  al.  (Ref.  27)  and  Cheret  and 
Barbier  (Ref.  28)  may  be  subject  to  a  different  interpretation. 

In  order  to  determine  the  possible  paths  for  alkali  excited  state  interactions,  we  have 
undertaken  a  series  of  ab  initio  calculations  of  the  electronic  states  of  Li2  up  to  the 

Li  [  ^  S(3s)]  +  Li  [  ^  S(3s)]  dissociation  limit.  Since  the  formation  of  the  ion-pair 

Li"*'  [  *  S]  +  Li"  [  *  S]  can  occur  only  for  ^2^  and  symmetries,  we  have  restricted  our 

studies  to  these  two  representations  (see  TSible  2). 

Theoretical  potential  curves  for  the  ^2^  and  symmetries  of  Li2  were  obtained  firom 

valence  configuration  interaction  (VCI)  calculations.  A  64  function  Slater  orbital  basis,  containing 
a  (2slp)  optimized  valence  basis  augmented  by  a  (3s3p3d2f)  Rydberg  basis,  was  used  in  this  study. 
The  basis  was  transformed  to  Dooh  symmetry  orbitals  and  Cl  calculations  were  performed  in  the 
space  of  all  symmetry  adapted  configurations  having  the  lOg  orbital  doubiy  occupied.  The 

resulting  Cl  expansions  contained  344  and  331  configurations  for  the  and  ^  symmetries, 

respective^.  As  the  interactions  in  this  system  are  very  long-range,  calculations  were  performed 
over  a  range  of  intemuclear  separations  from  5  to  50  bohrs. 

The  results  of  our  theoretical  calculations  are  shown  in  Fig.  5  and  Fig.  6  for  the  ^2^  and 

^2^  symmetries,  respectively.  The  ioiuc  curves  which  dissociate  to  Li"*"  [  ^  S]  +  Li"  [  ^  S(2s^)l 
are  clearly  evident  in  both  Fig.  5  and  Fig.  6,  where  they  exhibit  a  nearly  diabatic  crossing  behavior 
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with  the  normal  valence  states  of  Ia2  at  large  intemuclear  separations,  but  show  strong  mixing  at 
distances  ^  10  A .  Ion-pair  formation  may  occur  for  long-range  interactions  of  the  type 


U  [2S(2s)]  +  U**  [2S(n  =  5)]  -  U+  [^S]  +  U-[^S]  (28) 

but  clearly  there  are  few  effective  curve  crossing  channels  leading  to  Li'^  +  Li~ .  The  curve 
crossing  mechanism  suggested  by  Ciocca,  et  al.  (Ref.  27),  thus  appears  not  to  be  an  important 
mechanism  for  ion-pair  production  in  the  Li2  ^tem. 

A  second  mechanism  suggested  by  Ciocca  et  al.  (Ref.  27)  for  alkali  negative  ion  formation 
would  involve  excited  state  negative  ion  formation  via 

U*  (2p(2p)]  +  U**  (n  s  9)  -*  [*S]  +  U-[‘P(2s2p)]  (29) 

or  from  electron  attachment  to  a  Itydberg  state  of  li  *  *  followed  by  radiative  stabilization  to 
the  ground  Li~  [  ^  S(2s^)]  state.  Both  of  these  processes  can  be  ruled  out  for  the  Li2  system  since 
the  autoionization  lifetime  of  the  Li~  [^P(2s2p)]  state  is  known  to  be  very  short  (Ref.  89)  relative 
to  radiative  stabilization  and  higher  autoionizing  states  should  exhibit  even  shorter  lifetimes. 

As  mentioned  above,  McGeoch  and  Schlier  (^f.  29)  have  recently  found  that  the  most 
prominent  process  occurring  in  lithium  Itydberg  atom  collisions  is  associative  ionization: 

U  +  Li  Li2  +  e  (30) 

This  reaction  has  excess  exothermicity  beyond  the  dissociation  energy  of  Li^  ■  Some  mechanism 

for  disposing  this  excess  energy  must  be  developed,  either  by  radiation  or  by  internal  modes.  It  is 
very  unlikely  that  the  excess  energy  is  released  as  kinetic  energy  since  very  hot  electrons  would  be 
formed  which  would  result  in  secondary  ionization  processes.  An  alternate  possibility  is  that  an 
excited  metastable  state  of  Li^  formed  with  sufficient  lifetime  to  be  detected  and  mass 
analyzed. 

Finally,  a  careful  experimental  crossed  beam  study  of  He  Itydberg  atom-alkali  Itydberg 
atom  interactions  has  been  suggested  by  J.  Peterson  (SRI).  This  study  would  include  kinetic 
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energy  analysis  of  the  products  as  well  as  laser  selected  reactants.  The  main  reaction  appears  to 
form  a  direct  ion-pair  production  mechanism: 


He**  [^S(2s)]  +  K^ydbcrg  (Rb)  -  He’  +  K+(Rb-^) 


(31) 


The  interaction  potentials  and  energetics  of  such  excited  state  noble  gas-alkali  collisions  are 
unknown. 

7.  Dissociative  Attachment  Studies 

The  process  of  dissociative  electron  attachment  to  a  molecule  is  understood  (Ref.  90)  to 
proceed  via  the  intermediate  formation  of  an  electron-molecule  resonance  that  is  capable  of 
autodetaching  the  temporarily  bound  electron  with  a  finite  lifetime  that  is  related  to  the  width,  r , 
of  the  resonance.  The  potential  curves  of  the  neutral  molecule  (H2,  II2  or  LiH)  and  of  the 
resonant  state  ( He2>  Li2 ,  Li  H~ )  cross  at  an  intemuclear  separation  R  =  Rj  such  that,  for 
R  2  Rs,  the  autodetachment  of  the  electron  is  energetically  not  possible  and  the  resonance  turns 
into  a  stable  bound  state.  The  cross  section  for  dissociative  electron  attachment  is  determined  by 
the  asymptotic  value  of  the  radial  nuclear  wave  function  ^  (R)  of  the  resonant  state.  Using  energy 
normalized  continuum  functions,  we  have: 


Here  h  K/M  is  the  relative  veloci^  of  the  ion-atom  pair  after  the  attachment  process.  Note  that 
the  geometrical  cross  section,  ^/kf,  provides  an  upper  bound  (Ref.  91)  to  the  attachment  cross 

section  since  the  quantity  within  the  curly  parenthesis  can  be  interpreted  as  a  product  of  two 
probabilities,  one  for  the  capture  of  an  electron  to  form  the  resonant  state  and  the  other  (often 
called  the  survival  probability)  for  the  nuclei  in  the  resonant  state  to  separate  out  to  the  stabiliza¬ 
tion  radius  Rj.  Using  a  semiclassical  analysis  the  cross  section  for  dissociative  electron  attach¬ 
ment,  obA  >  can  also  be  written  as  (Ref.  92)  a  product  of  two  factors: 

obA  =  (33) 

The  first  factor,  Ocap ,  is  interpreted  as  the  cross  section  for  the  formation  of  a  resonant 
anion  state  by  the  capture  of  the  incident  electnxi  and  the  second  factor  S  (the  survival  probabili- 
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ty)  is  interpreted  as  the  probability  that  the  separation  of  the  nuclei  in  the  resonant  anion  state 
increases  up  to  the  stabilization  radius  Rj  without  electron  autodetachment  having  occurred.  This 
separation  ensures  that  the  process  of  dissociative  electron  attachment  will  occur.  The  survival 
probability  S,  in  the  semiclassical  analysis,  is  approximated  by 


exp(-rr/ft)  (34) 

where  r ,  the  time  taken  by  the  nuclei  to  separate  after  resonance  formation  to  the  stabilization 
R$,  is  directly  proportional,  by  simple  kinematical  considerations,  to  M^.  T  is  the  average  width 
of  the  resonance. 

The  computational  aspects  of  Eq.  (32)  consist  of:  a)  accurate  potential  energy  curves  for  the 
neutral  and  corresponding  negative  ion  species  over  a  range  of  intemuclear  separations  (R)  from 
the  repulsive  wall  to  separated  atom  states  and  b)  calculation  of  the  R  dependence  of  the  complex 
part  of  the  potential,  r  (R).  The  latter  calculation  is  relatively  difficult  and  usually  involves  some 
type  of  extrapolation  from  a  modified  nuclear  system  charge  or  a  direct  analysis  in  complex 
coordinates. 

The  UH  and  Li  H"  systems  were  studied  using  the  DIATOM  (Ref.  93)  molecular  structure 
code  with  a  [13a,  6 ;r,  26]  STO  basis.  These  STO  basis  calculations  were  required  to  examine 
the  relative  behavior  of  the  LiH  and  Li  H~  species  at  very  short  intemuclear  separations  and  to 
insure  correct  asymptotic  behavior.  The  non-linear  STO  parameters  were  optimized  to  describe 
the  asymptotic  negative  atomic  ions,  Li~  and  H~ .  The  variationally  determined  parameters  have 
been  previously  reported  (Ref.  25).  A  full  Cl,  with  frozen  core  for  U,  was  employed  to  insure 
correct  size  consistency  at  the  asymptote  limits.  The  lowest  two  molecular  states  of  Li  H~  are 
variationally  stable  at  large  intemuclear  separations.  As  the  intemuclear  separation  is  decreased, 
the  possibility  of  curve  crossing  with  the  neutral  ground  state  of  LiH  exists  and  the  anion  Cl 
framework  can  become  variationally  unstable,  lb  treat  this  problem,  we  empltty  a  modified 
nuclear  charge  Hamiltonian,  similar  to  that  described  by  Nestmann  and  Peyerimhoff  (Ref.  94). 

The  modified  Hamiltonian  causes  the  anion  states  to  lower  below  the  corresponding  neutral 
molecule.  A  systematic  extrapolation  of  the  modified  Hamiltonian  is  then  employed  to  predict  the 
location  of  the  molecular  anion  states  in  the  autodetaching  regions. 

In  order  to  examine  the  stability  of  Li^Hy  clusters,  ab  initio  calculations  were  carried  out  for 
several  species  at  the  MP2  level  of  theoiy.  The  basis  set  chosen  was  the  Gaussian  6-311G  triple 
split-valence  set  (Ref.  95),  augmented  by  d-polarization  functions  [as 0.20]  for  Li  and  p-polariza- 
tion  [a  s  0.75]  for  H.  In  addition,  a  diffuse  s-function  [a =0.036]  and  a  diffuse  sp-shell 
[a =0.0074]  were  added  to  hydrogen  and  lithium  respectively,  to  better  describe  the  negative  ion 
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charge  distributions.  All  calculations  were  performed  using  CADPAC  (Ref.  96),  an  electronic 
structure  code  that  can  perform  geometry  optimization  and  calculate  force  constant  matrices 
analytically.  Optimized  geometries  were  calculated  at  the  MP2  level  of  theory.  Harmonic  vibra¬ 
tional  frequencies  were  subsequently  calculated  at  the  MP2  stationary  points. 

The  results  of  our  calculations  for  the  LiH  and  LiH~  systems,  using  the  DIATOM  code,  are 
shown  in  Fig.  7  for  states  of  symmetry  of  Li  H” .  In  Fig.  8,  we  illustrate  the  short  range 

behavior  of  the  lowest  and  first  excited  ^2"^  states  of  Li  H" .  It  is  evident  that  the  ground 

state  of  Li  H~  is  bound,  relative  to  of  LiH,  for  all  intemuclear  separations.  Our  calculated 

adiabatic  electron  affinity  is  0.33  eV,  in  good  agreement  with  previous  theoretical  studies  (Ref.  97). 
The  first  excited  ^  S'*’  state  of  Li  H' ,  which  asymmetrically  correlates  to  Li" + H,  exhibits 

repulsive  behavior  in  the  region  3.0  R  ^  6.0  A .  Hius,  dissociative  attachment  of  e  -I-  LiH 
mainly  forms  Li"  for  Ecoii  2:  3.0  eV  The  formation  of  Li  -H  H“  products  may  occur  by  non-adia- 

batic  coupling  of  the  X^S^  state  of  Li  H"  to  the  continuum  of  e  -i-  LiH  [  X^S"*^  ]  for  Ecoii  ^  2.1 
eV  However,  the  cross  section  for  such  an  indirect  electron  capture  mechanism  is  predicted  to  be 
very  small. 

Figure  9  shows  a  comparison  of  the  potential  curves  of  the  ground  states  of  the  molecules 
H2,  LiH  and  Li2  as  well  as  the  potential  curves  of  the  states  of  the  corresponding  molecular 
anions  that  are  primarily  responsible  for  dissociative  electron  attachment  to  these  molecules.  As 
one  progresses  from  H2  to  LiH  to  Li2,  the  equilibrium  intemuclear  separation  increases  from  1.40 
a.u.  to  3.08  a.u.  to  5.06  a.u.  and  the  dissociation  energy  of  the  molecule  decreases  from  4.73  eV  to 
2.33  eV  to  1.06  eV.  Furthermore,  the  energy  spacing  between  the  same  two  low-lying  adjacent 
vibrational  levels  decreases  as  one  proceeds  from  H2  to  UH  to  Li2,  in  such  a  manner  that  the 
total  number  of  bound  vibrational  levels  increases  on  increasing  the  molecular  mass.  These 
energetics  have  different  quantitative  effects  on  the  threshold  for  dissociative  electron  attachment 
to  these  molecules  as  well  as  on  the  factor  by  which  the  attachment  cross  section  is  enhanced  on 
vibrationally  exciting  these  molecules.  The  threshold  energy  for  dissociative  attachment  to  a 
molecule  that  is  rovibrationally  excited  to  a  particular  level  (v4)  is  Eth^^  =  D-EA-Eyj,  if  Eyj 
<  D-EA,  and  Etj,®'^ =0  otherwise.  The  lower  energy  threshold  for  dissociative  attachment  thus 
enhances  the  attachment  cross  sections  as  we  proceed  from  H2  to  LiH  to  Li2. 

The  possible  role  of  Li^Hy  clusters  was  examined  through  ab  initio  calculations  of  Li2H, 

Li^,  Li2H2  and  their  respective  negative  ion  states.  All  clusters  calculations  were  carried  out 
using  the  CADPAC  code.  Optimized  geometries  were  calculated  at  the  MP2  level  of  theory. 
Harmonic  vibrational  analyses  were  subsequently  computed  at  the  MP2  stationaiy  points.  In 
several  instances,  the  cluster  anion  is  unbound  relative  to  the  corresponding  neutral  molecule.  In 
such  cases,  a  check  on  vibrational  stability  is  required  to  insure  that  there  are  no  instabilities 
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toward  anion  decomposition.  Further,  the  endothermicity  of  the  anion  calculation  must  be  studied 
as  a  function  of  the  addition  of  diffuse  basis  functions.  The  anions  studied  here  did  not  exhibit  a 
strong  energy  dependence  on  the  diffuse  functions,  suggesting  that  they  may  represent  Feshbach 
structures.  This  feature,  however,  was  not  explored  in  detail  and  some  uncertainty  can  be  assigned 
to  our  energetics  of  the  thermodynamically  unstable  anion  clusters.  In  Fig.  10,  we  illustrate  the 
stationary  geometry  of  Li2H  and  the  corresponding  negative  ion.  We  find  that  the  C2v  structure 
for  Li2H  has  vibrational  stability  and  is  thermodynamically  stable  relative  to  Li2  +  H2.  In 
addition,  the  corresponding  C2v  anion  is  slightly  bound  (0.06  eV)  relative  to  the  neutral  species. 
These  results  are  in  qualitative  agreement  with  previous  theoretical  studies  of  Li2H  by  Cardelino, 
et  al.  (Ref.  98),  that  were  carried  out  at  the  SCF  level  of  theory. 

The  Li3H  species  was  examined  both  as  a  C3V  (trigonal  pyramid)  and  as  a  C2v  (kite-like) 
structure.  Our  previous  studies  (Ref.  99)  of  Li3H  indicated  that,  although  both  structures 
exhibited  stable  vibrational  frequencies,  only  the  C2v  structure  was  thermodynamically  stable.  In 
Fig.  11,  we  illustrate  the  two  stable  geometries  of  neutral  Li3H.  Of  these,  the  C2v  structure 
represents  the  global  minimum  for  the  system.  The  Li3  H~  anion  is  not  stable  as  a  C3V  structure 
and  collapses  to  Li3  +  H"  .  The  C2v  structure  of  Li3  H”  has  the  stationary  structure  shown  in 
Fig.  11,  but  a  vibrational  analysis  indicates  that  this  is  probably  a  saddle  region  on  the  potential 
energy  surface.  In  addition,  this  region  is  thermodynamically  endothermic  by  0.24  eV,  relative  to 
the  neutral  Li3H. 

Finally,  calculations  were  carried  out  for  the  (LiH)2  dimer,  as  a  linear  (Coov)>  rhomboid 
(D2h)  and  y-structure  (C2v)-  The  results  are  shown  in  Fig.  12  where  we  show  that  the  linear 
conformation  has  one  imaginary  frequency,  and  thus  represents  a  saddle  point  on  the  surface. 

The  C2v  structure  is  totally  unstable  and  collapses  to  separate  Li2  +  H2  molecules.  The  only 

vibrationally  stable  geometry  is  the  ^  Ai  state  as  a  D2h  structure.  The  corresponding  negative 

ion  Li2  H2  is  vibrationally  unstable  and  unbound  by  over  2.0  eV,  relative  to  Li2H2(D2h). 

A  summary  of  the  thermodynamics  of  these  UxHy  clusters  is  given  in  Table  3.  All  energies 
are  given  relative  to  Li2  tmd  H2  as  gas  phase  species.  As  previously  known,  the  formation  of  LiH 
is  endothermic  by  0.740  eV.  However,  higher  clusters  are  all  thermodynamically  stable;  Li2H 
(-0.280  eV),  Li3H  (-1.034  eV)  and  y2H2  (-1263  eV). 

The  trend  of  our  studies  to  date  indicates  that  large  LixHy  (and  by  analogy,  CsxHy)  clusters 
are  thermodynamically  stable.  In  particular,  the  Li2H2  species,  as  a  C2v  structure,  may  be  an 
important  component  of  alkali-hydrogen  mixtures.  This  sp)ecies  can  dissociatively  attach  an 
electron  to  form  Li2H  +  H"  for  Ecoii  ^  2.0  eV  The  Li3H  (C2v)  species  should  also  exhibit  DA  to 
form  Li3  +  H~ ,  but  the  concentration  of  this  molecule  will  be  lower  than  that  of  the  more  stable 
Li2H2  cluster.  The  mechanism  for  DA  to  these  LixHy  clusters  can  be  understood  by  following 
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their  intrinsic  reaction  coordinate  pathways  leading  to  dissociation.  In  addition,  the  role  of  higher 
order  clusters  in  volume  DA  processes  is  uncertain. 

8.  Studies  of  Negative  Ion  Thermochemistry 

In  recent  years,  the  dissociation  energy  of  a  C-H  bond  of  acetylene  HCCH-»H  +  CCH  has 
been  the  subject  of  a  number  of  studies,  both  experimental  (Refs.  100-102)  and  theoretical 
(Ref.  103).  All  of  these  studies  were  in  close  agreement,  reporting  values  or  upper  bounds  for 
Do(H-CCH)  ranging  from  131.9  kcal/mol  (Ref.  100)  to  133.5  kcal/mol  (Ref.  103).  However,  two 
very  convincing  new  measurements  (Refs.  104  and  105)  of  this  dissociation  energy  Dq  <  126.647 
±  0.002  kcal/mol  (Ref.  104)  and  Dq  127  ±  1.5  kcal/mol  (Ref.  105)  are  in  close  agreement  with 
Benson’s  thermochemical  estimate  (Dq  126  ±  1  kcal/mol  (Ref.  106))  and  have  raised  doubts 
about  the  accuracy  of  the  higher  values.  In  an  attempt  to  determine  which  value  is  correct,  we 
have  performed  a  new  ab  initio  calculation  of  this  dissociation  energy  including  extrapolation  to 
the  complete  basis  set  (CBS)  limit  (Refs.  107-115).  Our  analysis  indicates  that  the  probable  error 
in  the  new  calculation  is  less  than  1  kcal/mol.  Our  calculated  dissociation  energy  (Dq  131.54 
kcal/mol)  is  consistent  with  the  higher  values,  but  not  with  the  most  recent  experiments  (Refs.  104 
and  105). 

In  conjunction  with  this  analysis  of  the  thermochemistry  of  C2H2,  we  have  calculated  the 
electron  affinity  of  the  C2H  radical.  The  recent  accurate  measurement  of  the  electron  affinity  of 
CCH  (2.969  ±  0.010  eV)  by  Ervin  et  al.  (Ref.  116)  offers  an  additional  opportunity  to  assess  the 
accuracy  of  our  calculated  results.  If  there  is  a  systematic  error  in  our  calculated  CCH  energy,  it 
would  also  affect  our  prediction  of  the  CCH  electron  affinity.  Using  the  basis  set  described  above, 
augmented  by  an  additional  diffuse  carbon  d  function,  we  find  EA(CCH)  =  2.961  ±  0.015  eV,  in 
excellent  agreement  with  Ervin  et  al. 

Subsequently,  Ervin  and  Lineberger  (Ref.  117)  reported  new  experimental  values  for  the 
electron  affinity  (2.%9  ±  0.006  eV)  and  dissociation  energy  (187.2  ±  2.5  kcal/mol)  of  C2~ 

(Ref.  117).  As  a  further  test  of  our  methodology,  we  have  performed  CBS-QCI  calculations  for 
comparison  with  these  experimental  results.  The  calculated  results  are  EA  (C2)  =  2.92  ±  0.02  eV 
Do  (C2“)  =  188.9  ±  1.5  kcal/mol  and  Do  (C2)  =  143.7  ±  1.5  kcal/mol.  The  calculated  values  for 
EA  (C2)  and  Do  (C2~)  are  found  to  be  in  good  agreement  with  the  experimental  results,  although 
the  errors  slightly  exceed  our  expected  values.  One  interesting  point  emerges  from  further 
comparison.  Using  a  thermodynamic  cycle,  Ervin  and  Lineberger  derive  Do  (C2~)  from  their 
measured  value  of  EA  (C2)  and  the  accurately  known  value  of  EA  (C),  and  Do  (C2).  For  Do  (C2), 
which  their  experiment  does  not  determine,  they  use  a  value  of  141.0  ±  2.5  kcal/mol,  taken  from 
literature  values.  A  higher  value  of  Do  (C2),  such  as  the  143.2  kcal/mol  recommended  by  Herzberg 
and  Huber  (Ref.  118)  would  be  in  better  agreement  with  our  calculated  Do  (C2),  and  in  combina¬ 
tion  with  the  accurate  experimental  electron  affinities,  would  give  a  value  (189.4  kcal/mol)  for 
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Do  (C2  )  that  is  also  in  better  agreement  with  our  theoretical  value.  Thus,  these  calculations 
suggest  that  the  accepted  value  of  Do  (C2)  is  roughly  2-3  kcal/mol  too  low. 

CBS-QCI  studies  have  also  been  performed  on  the  vinylidene-acetylene  isomerization 
process,  which  has  been  recently  observed  by  Lineberger  and  co-workers  (Ref.  119).  The  calcu¬ 
lated  isomerization  barrier  (AEe  =  2.2  ±  0.5  kcal/mol)  and  energy  change  from  vinylidene  to 
acetylene  (AEo  =  -43.9  ±  0.5  kcal/mol)  are  in  excellent  agreement  with  recent  experimental 
values  (AEe  =  2  kcal/mol,  and  AEo  =  -43  ±  2  kcal/mol).  In  spite  of  the  small  barrier  height 

for  this  strongly  exothermic  reaction,  the  transition  state  is  located  halfway  between  the  reactant 
and  product  in  apparent  violation  of  the  Hammond  postulate  (Ref.  120).  We  have  shown  tb  ;ne 
potential  energy  surface  for  the  isomerization  can  best  be  understood  by  considering  two  distinct 
processes.  The  location  of  the  transition  state  is  determined  by  the  hydrogen  migration  which  has 
a  significant  barrier,  whereas  the  exothermicity  results  from  the  conversion  of  the  lone  pair  of 
vinylidene  to  a  ;r  bond  in  acetylene.  Each  of  these  processes  individually  satisfies  the  Hammond 
postulate. 
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TABLE  2 

LOW-LYING  MOLECULAR  STATES  OF  Li.  AND  THEIR  DISSOCIATION  LIMITS 


Separated  atoms  E(cm~l) 


Li  Li 

2S  (2s)  +  2s  (28)  0.0 

o  o 

2Sg(28)  +  2p^(2p)  14903.0 

2S  (2s)  +  2s Os)  27206.0 

o  o 

2p^(2p)  +  2p^(2p)  29807.0 


2Sg(28)  +  2p^(3p)  30925.0 


2S  (28)  +  2d  (3,1)  31283.0 

o  o 


LiJl^J^]  +  e 

33200,0 

2Sg(28)  +  2s^(48) 

35012.0 

2Sg(28)  +  2p^(4p) 

36469.0 

C.  E.  Moore,  Atomic  Energy  Levels, 


Molecular  states 

iz+d),  3z;;;(i) 

Ir^d).  aijd),  3e;^(i) 

iRgd),  in^d),  ^Hgd), 

iz+d),  lE+d),  3eJ(i).  3z;^(i) 
1e+(2),  3z+(2),  ^Egd) 

iRgd),  3n^d),  ^ng(i),  ^n^d) 
lAgd),  X(i), 
lEjd),  lE^d).  3j:+(i), 

3z^d), 

3ng(i),  3ny(i) 

izjd),  lE+d),  3z+(l), 

3e^(i),  ^Hgd),  ^n^d),  ^Hgd), 

3n^(l),  lAgd),  lA^d),  3Ag(l), 

lE+d).  3e+(1),  3j;+(i) 
lE+d),  1e;J(1),  3j;+(1), 

3e+(i).  iRgd),  in^^d), 

3ng(i),  3n^^(i) 

(1949). 
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TABLE  2  (Continued) 

2Sg(2s)  +  2Dg(4d) 

36623.0 

'ijd), 

‘ijd). 

3z*d). 

‘Bgd). 

‘"ud), 

V”- 

■  X“>- 

‘V*>- 

‘4„d), 

V‘>> 

2Sg(28)  +  2F^(4f) 

36630.0 

HIM. 

‘Ejd). 

HIM. 

‘V'>- 

iii„d), 

>i,d), 

‘V”> 

3A^d), 

X(». 

H^M. 

IVI). 

3.gd), 

2S  (2s)  ♦  2s  (5s) 

O  O 

38299.0 

HIM. 

HIM. 

3!:*d). 

3E;d) 

Li+  [Is  (ls2)]  + 

Li-  I1S*(2s2)] 

38474.0 

him. 

H^M 

2Sg(28)  ♦  2p^(5p) 

39015.0 

HIM. 

>V‘>. 

3l*d), 

■n,d). 

3n„d), 

V‘>- 

3V 1) 

2Sg(2s)  +  2D^(5d) 

39094.0 

him. 

■i^d). 

3z*d). 

H*M. 

■n^d), 

>ll„d). 

3llgd), 

3ll„(l). 

\M. 

‘4„d). 

3Agd). 

X”) 

2Sg(2s)  +  2F^(5f) 

39104.0 

'ijd). 

h*m. 

3l*d), 

3V1). 

‘V”> 

in^(i), 

%»>• 

3n„d), 

>V‘). 

lA^(l). 

V‘>- 

3A„d). 

l#u(l). 

3^1). 

3«„d) 

2Sg(2s)  ♦  Doublets 
(5<n<9) 

39987.0  ■*. 

42003.0 

singlet , 

E  9  n  9  A  9 

triplet 

(g.u) 
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2p^(2p)  +  2Sg(3s) 


2s  (28)  +  Doublet 
(9<n<«) 

Li+  [Is  (ls2)]  + 
Li  l2Sgr28)] 

2p^(2p)  +  2p^(3p) 


2p^(2p)  +  2D^(3d) 


2p^(2p)  +  Doublet 
(n*4,  5) 

2S  (38)  +  2s  (3s) 

O  O 

Li-^  [Is  (182)]  + 
Li'  [lp*(282p)l 


TABLE  2  (Concluded) 

42110.0 

3e+(1),  ifigd),  X(l). 

3n^(i),  3n^(i) 

42118.0  Singlet,  triplet 

Z,  n,  A  (g,u) 

43482.0  2j:+(i),  2j;+(i) 

o 

45829.0  ^^g^2),  1e+(2).  3j;+(2), 

3e*(2),  ^Egd).  ^J'd), 

3e'(1),  3j;'(1),  ^ng(2), 

llly(2),  3ng(2),  3n^(2),  ^Ag(l), 
lA^jd),  ^Agd),  3A^d) 

46187.0  lE*(2),  ^Zg(l). 

IZ'd).  3e+(2),  3j:+(2). 

3Eg(l),  3j:-(1),  lng(3), 

llly(3),  3llg(3),  3n^(3),  lAg(2), 

lAy(2),  ^Ag(2),  3/^(2)  ^♦gd). 

^•yd),  3j^(l),  3j^(l) 

51088.0  -►  Singlet,  triplet 

53782.0  E,  n.  A,  ♦,  T  (g,u) 

54412.0 

57920.0  ^j:*d),  ^lyd),  in_(i), 

In^d) 
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TABLES 

THERMODYNAMICS  OF  U*Hy  CLUSTERS 


Reaction 

AE(eV) 

Li2(g)  +  H2(g)^2UH(g) 

+0.740 

(1) 

U2(g)+  l/2H2(g)-.U2H(C2v) 

-0.280 

(2) 

3/2  U2  (g)  +  1/2  H2  (g)  ^  y  3H  (C2v) 

-1.034 

(3) 

Li2(g)  +  H2(g)-.U2H2(D2h) 

-1.263 

(4) 
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FIG.  7 
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FIG.  8 
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FIG.  11 
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APPENDIX  I-l 


POTENTIAL  ENERGY  CURVES  FOR  DISSOCIATIVE 
RECOMBINATION  OF  HeH+ 

H.H.  Michels 

United  Technologies  Research  Center 
East  Hartford,  Connecticut  06108 


ABSTRACT 

The  potential  energy  curves  pertinent  to  the  process  of  dissociative  recombination  (DR)  of  e  + 
HeH'*’  have  been  calculated  using  an  ab  initio  valence  configuration  interaction  (VCI)  method. 
Recombination  into  both  the  singlet  and  triplet  ion  channels  has  been  considered  through  analysis  of 
the  low-lying  doublet  and  quartet  states  of  HeH.  Wth  the  exception  of  the  ground  state,  all  of  the 
doublet  states  up  to  the  HeH"*"  limit  are  known  to  be  Rydberg  in  character.  Similarly,  we  find  weak 

binding  for  the  lowest  quartet  states  of  HeH,  which  appear  to  be  mainly  Rydberg  to  the  a^Z^  state  of 
HeH'''.  The  nature  of  a  second  quartet  series,  which  should  be  Rydberg  to  the  lowest  ^11  and  first 

excited  ^Z'''  states  of  HeH''",  is  uncertain,  owing  to  possible  strong  perturbations  with  the  ion-pair 
states  of  HeH. 


Supported  in  part  by  AFOSR  under  Contracts  F49620-85-C-(K)95  and  F49620-89-C-0019. 
Published  as  a  chapter  in  Dissociative  Recombination:  Theory,  Experiment  and  Application,  eds. 
J.  B.  A.  Mitchell  and  S.  L  Guberman,  World  Scientific  Pub.  Co.  Inc.,  1989,  pp.  97-108. 
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APPENDIX  1-2 

THEORETICAL  STUDY  OF  IONIZATION  AND  ELECTRON 
ATTACHMENT  PROCESSES  IN  RYDBERG  ATOM 
COLLISIONS 


H.H.  Michels  and  R.  H.  Hobbs 
United  Technologies  Research  Center 
East  Hartford,  Connecticut  06108 


ABSTRACT 

Theoretical  studies  of  collisional  ionization  and  electron  attachment  mechanisms  involving  light 
elements  (H,  He,  Li)  are  being  carried  out  in  support  of  particle  beam  sources  and  beam  neutralization 
concepts.  In  particular,  the  electronic  structure  of  stable  anions,  their  resonance  and  metastable 
states,  and  radiative  and  collisional  processes  describing  their  formation  and/or  destruction  are  being 
investigated.  The  volume  production  of  atomic  anions  via  electron-molecule  dissociative  attachment 
has  previously  been  examined,  both  experimentally  and  theoretically,  for  e  +  Hz  and  e  +  Liz 
collisions.  An  alternate  production  route,  involving  ion-pair  formation  from  Rydberg-Rydberg  atom 
collisions,  is  currently  being  investigated.  The  possible  mechanisms  include:  ionic-Rydberg  molecule 
curve  crossings,  collisional  radiative  stabilization,  and  associative  ionization  followed  by  secondary 
electron-neutral  attachment  processes.  Recent  experiments  by  McGeoch,  et  al.  indicate  that 
associative  ionization  is  highly  favored  over  direct  ion-pair  production  in  the  Li  system.  However, 
both  Ciocca,  et  al  and  Cheret  and  Barbier  report  efficient  ion-pair  production  in  collisions  between 
heavier  (Na,  Rb)  alkali  Itydberg  atoms.  Several  p(K;sible  reaction  mechanisms  have  been  examined  for 
these  systems. 


Supported  in  part  by  AFOSR  under  Contracts  F4%20-85-C-0095  and  F49620-89-C-0019. 
Published  in  SPIE-^crowave  and  Particle  Beam  Sources  and  Directed  Energy  Concepts,  Vol. 
1061, 1989,  pp.  597-604. 
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APPENDIX  1-3 

LASER  INDUCED  CHEMILUMINESCENCE  OF  THE  LiMg 

EXCIMER 

G.  Pichler,  A.  M.  Lyyra,  P.  D,  Kleiber, 

W.  C.  Stwalley,  R.  Hammer,  K.  M.  Sando 
Center  for  Laser  Science  and  Engineering 
University  of  Iowa 
Iowa  City,  lA  52242  USA 

H.  H.  Michels 

United  Ibchnologies  Research  Center 
East  Hartford,  CT  06108 


ABSTRACT 

A  UV  Ar^  laser  was  used  to  excite  specific  levels  in  the  C^Ilu  state  of  the  ^  Li2  dimer.  A 
subsequent  photochemical  reaction  with  atomic  Mg  produced  the  LiMg  excimer  in  excited 
C^n  and  states.  The  observed  bound-free  »ccimer  emission  spectrum  was  compared  with 
theoretical  emission  profiles  based  on  potential  energy  curves  and  transition  dipole  moments 
obtained  from  ab  initio  calculation.  These  calculations  also  indicate  avoided  crossings  between  the 
A^n  and  C^H  states,  and  and  D^Z'*’  states  of  LiMg. 


Supported  in  part  by  AFOSR  under  Contracts  F49620-85-C-0095  and  F4%20-89-C-0019. 
Published  in  Chemical  Physics  Letters,  156, 1989,  pp.  467-471. 
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APPENDIX  1-4 

GAS  PHASE  ALKALI-HYDROGEN  INTERACTIONS 
IN  NEGATIVE  ION  SOURCES 

H.H.  Michels 

United  Technologies  Research  Center 
East  Hartford,  Connecticut  06108 

J.  M.  \^dehra 
>\^yne  State  University 
Detroit,  MI  48202 


ABSTRACT 

The  role  of  LixHy(CsxHy)  molecules,  which  can  be  formed  from  seeding  an  alkali  into  a  hydrogen 

plasma,  is  examined  through  ob  initio  calculations  of  the  structure  and  stability  of  such  species.  The 
simplest  alkali  hydride,  UH,  supports  a  bound  anion  for  all  intemuclear  separations.  Larger 
structures,  such  as  Li2H2  and  UbH  ,  do  not  support  stable  anions  but  exhibit  thermodynamic 

stability  as  gas  phase  molecules.  Their  possible  roles  in  dissociative  attachment  to  form  H~  and/or 
Li~  is  examined. 


Supported  in  part  by  AFOSR  under  Contract  F4%20-89-C-0019. 

Published  in  Particle  and  Fields  Series  40,  Proceedings  of  the  Fifth  International  Symposium  on 
the  Production  and  Neutralization  of  Negative  Ions  and  Beams,  Conf.  Proc.  No.  210  (Brookhaven, 
NY,  1989),  ed.  Ady  Hershcovitch,  p.  142, 1990. 


1-5 


R92-928101 


APPENDIX  1-5 

DISSOCIATIVE  ELECTRON  ATTACHMENT  TO  LIGHT  MOLECULES:  A 
COMPARATIVE  STUDY  OF  H2,  LiH  AND  U2 

H.H.  Michels 

United  Technologies  Research  Center 
East  Hartford,  Connecticut  06108 

J.  M.  W^dehra 

Department  of  Physics  and  Astronomy 
Wayne  State  University 
Detroit,  MI  48202 

ABSTRACT 

In  this  paper  we  have  compared  the  energetics  of  three  light  molecules  (H2,  LiH  and  Li2)  to 
form  negative  ions  by  the  process  of  dissociative  electron  attachment.  For  two  light  homonuclear 
molecules  (H2  and  Li2)  we  have  done  explicit  calculations,  using  the  local  formalism  of  the  resonance 
model,  to  obtain  the  electron  attachment  cross  sections  and  the  rates  of  negative  ion  formation. 


Supported  in  part  by  AFOSR  under  Contract  F4962ft-89-C-0019. 

Published  in  Particles  and  Fields  Series  40,  Proceedings  of  the  Fifth  International  Symposium  on 
the  Production  and  Neutralization  of  Negative  Ions  and  Beams,  Conf.  Proc.  No.  210  (Brookhaven, 
NY,  1989),  ed.  Ady  Hershcovitch,  p.  114, 1990. 


1-6 


R92-928101 


APPENDIX  1-6 


ON  THE  C-H  BOND  DISSOCIATION  ENERGY  OF  ACETYLENE 

J.  A  Montgomery,  Jr. 

United  Technologies  Research  Center 
East  Hartford,  Connecticut  06108 


G.  A  Petersson 

Hall-Atwater  Laboratories  of  Chemistry 
Wesleyan  University 
Middletown,  CT  06457 


ABSTRACT 

Ab  initio  calculations,  including  extrapolation  to  the  complete  basis  set  limit,  are  reported  for  the 
C-H  bond  dissociation  energy  of  acetylene.  These  calculations,  which  have  a  RMS  error  of  0.51 
kcal/mol  per  bond  for  the  atomization  energies  of  13  well-known  molecules,  give  DofH-CCH)  * 
131.54  kcal/mol,  in  good  agreement  with  previous  results,  but  about  5  kcal/mol  above  the  recent 
measurement  of  Green,  Kinesey,  and  Field  (Dq  ^  126.647  ±  0.002  kcal/mol),  and  that  of  Segall,  Lavi, 
Wen,  and  Wittig  (Dq  =  127  ±  1.5  kcal/mol).  The  electron  affinity  of  CCH  is  calculated  to  be  within 
0.18  kcal/mol  of  the  recent  measurement  of  Ervin  et  al.  (2.969  ±  0.010  eV). 


Supported  in  part  by  AFOSR  under  Contract  F49620-89-C-0019. 
Published  in  Chemical  Physics  Letters,  1^  1990,  pp.  75-78. 
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APPENDIX  1-7 

STRUCTURE  AND  STABILITY  OF  Li^H,  MOLECULES 

AND  ANIONS 


H.H.  Michels 

United  Technologies  Research  Center 
East  Hartford,  Connecticut  06108 

J.  M.  Wadehra 
W^yne  State  University 
Detroit,  MI  48202 


ABSTRACT 

The  role  of  LixHy  (CsxHy)  molecules,  which  can  be  formed  from  seeding  an  alkali  into  a  hydrogen 
plasma,  is  examined  through  ab  initio  calculations  of  the  structure  and  stability  of  such  species.  The 
simplest  alkali  hydride,  LiH,  supports  a  bound  anion  for  all  intemuclear  separations.  Larger 
structures,  such  as  Li2H2  and  Li3H,  do  not  support  stable  anions  but  exhibit  thermodynamic  stability 
as  gas  phase  molecules.  Their  possible  roles  in  dissociative  attachment  to  form  H~  and/or  Li~  is 
examined. 


Supfwrted  in  part  by  AFOSR  under  Contract  F49620-8^-C-0019. 

Published  in  Ae  International  Journal  of  Quantum  Chemistry,  Sanibel  Symposia  Supplement, 
1990,  pp.  521-530. 
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APPENDIX  1-8 

COMMENT  ON  “AN  EXPERIMENTAL  AND 
THEORETICAL  STUDY  OF  THE  NEGATIVELY  CHARGED 
HELIUM  DIMER  He^ 

H.  H.  Michels,  R.  H.  Hobbs 
United  Technologies  Research  Center 
East  Hartford,  CT  06108  USA 

Y.  K.  Bae  and  J.  R.  Peterson 
SRI  International 
Menlo  Park,  CA  94025  USA 

ABSTRACT 

It  is  shown  that  the  study  of  the  negatively  charged  helium  dimer  reported  by  Krishnamurthy 
and  Mathur  contains  serious  errors  in  both  the  theoretical  interpretation  and  experimental 
analysis  of  this  ion. 


Supported  in  part  by  AFOSR  under  Contract  F49620-89-C-0019. 
Published  in  Chemical  Physics  Letters,  170, 1990,  pp,  411-412. 
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APPENDIX  1-9 

A  COMPLETE  BASIS  SET  MODEL  CHEMISTRY.  III.  THE 
CBS-QCI  FAMILY  OF  METHODS 

G.  A.  Petersson  and  Thomas  G.  Tensfeldt 
Hall- Atwater  Laboratories  of  Chemistry 
Wesleyan  University 
Middletown,  CT  06457 

J.  A.  Montgomery,  Jr. 

United  Technologies  Research  Center 
East  Hartford,  Connecticut  06108 


ABSTRACT 

The  major  source  of  error  in  most  ab  initio  calculations  of  molecular  energies  is  the  truncation  of 
the  one-electron  basis  set.  A  family  of  complete  basis  set  (CBS)  quadratic  Cl  (QCI)  model  chemistries 
is  defined  to  include  corrections  for  basis  set  truncations  errors.  These  models  are  basis  sets  ranging 
from  the  small  6-31GttDZ+ P  size  basis  set  to  the  very  large  (14s9p4d2f,6s3pld)/[6s6p3d2f,4s2pld] 
atomic  pair  natural  orbital  basis  set.  When  the  calculated  energies  are  compared  with  the 
experimental  energies  of  the  first-row  atoms  and  ions  and  the  first  row-diatomics  and  hydrides:  H2, 
LiH,  Li2.  CH4,  NH3,  H2O,  HF,  liF,  N2,  CO,  NO,  O2,  and  F2,  two  very  promising  new  model  chemistries 
emerge.  The  first  is  of  comparable  accuracy,  but  more  than  ten  times  the  speed  of  the  G1  model  of 
Pople  and  co-workers.  The  second  is  less  than  one  tenth  the  speed  of  the  G1  model,  but  reduces  the 
RMS  errors  in  IP’s,  EA’s,  and  Dq’s  to  0.033  eV,  O.OD  eV  and  0.53  kcal/mol  per  bond  respectively. 


Supported  in  part  by  AFOSR  under  Contract  F4%20-89-C-0019. 
Published  in  Journal  of  Chemical  Physics,  94, 1991,  pp.  6091-6101. 
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APPENDIX  1-10 

TEMPERATURE,  KINETIC  ENERGY,  AND  ROTATIONAL 
TEMPERATURE  EFFECTS  IN  REACTIONS  INVOLVING  ISOTOPES 


A.  A,  Mggiano,  R.  A.  Morris,  J.  M.  Van  Doren  and  J.  F.  Paulson 
Phillips  Laboratory,  Geophysics  Directorate 
Ionospheric  Effects  Division  (LID) 

Hanscom  AFB,  MA  01731-5000 

H.  H.  Michels  and  R.  H.  Hobbs 
United  Ibchnologies  Research  Center 
East  Hartford,  Connecticut  06108 

ABSTRACT 

Data  on  four  reactions  involving  isotopes  taken  in  a  variable  temperature  selected  ion  flow  drift 
tube  are  presented.  A  study  of  the  reaction  of  0~  with  N2O  indicates  that  the  reaction  proceeds 
preferentially  by  bonding  of  the  0~  to  the  central  nitrogen  in  N2O.  The  preference  for  O'  attack  at  the 
central  nitrogen  over  attack  at  the  terminal  nitrogen  decreases  at  higher  temperatures.  In  the  atom 
abstraction  reaction  of  O'**  with  HD,  OH'*'  is  formed  more  efficiently  than  is  OD'*'  at  low 
temperatures  and  moderate  energy.  The  branching  fraction  favoring  OH  production  is  also  sensitive 
to  the  rotational  temperature  of  the  HD.  The  results  for  this  reaction  are  consistent  with  a  model  of  the 
reaction  based  on  the  long  range  part  of  the  ion-neutral  potential.  Rate  constants  for  the  reactions  of 
O'  with  H2,  D2,  and  HD  vary  with  mass  of  the  hydrogen  molecule  as  predicted  from  the  collision  rate 
constant  dependence.  These  reactions  proceed  by  two  channels:  hydrogen  abstraction  and  associative 
electron  detachment.  The  rate  constant  for  the  minor  hydrogen  abstraction  channel  increases  with 
increasing  kinetic  energy.  The  efficiency  of  the  abstraction  reaction  of  O'  with  D2  is  significantly 
smaller  than  that  in  the  reaction  with  H2.  More  OD~  as  compared  with  OH'  is  produced  in  the  reaction 
of  HD.  The  results  are  explained  by  a  direct  two  step  mechanism.  A  large  preference  for  OH~  over 
OD'  formation  in  the  atom  abstraction  reaction  of  O' with  CH2D2  is  observed.  An  ab  initio  potential 
for  the  reaction  path  reveals  a  barrier  and  indicates  that  zero  point  energy  effects  play  a  major  role  in 
the  observed  isotope  effects.  The  four  reactions  studied  in  this  paper,  when  taken  together,  show  that 
isotopically  labelled  reactants  can  play  a  major  role  in  elucidating  widely  different  reaction 
mechanisms. 


Supported  in  part  1^  AFOSR  under  Contract  F49620-89-C-0019. 
Submitted  for  publication  to  ACS  volume  on  Isotopic  Effects  in  Chemical 
Reactions  and  Photodissociation  Processes,  1991. 
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APPENDIX  I-ll 

VINYLIDENE  AND  THE  HAMMOND  POSTULATE 

G.  A.  Petersson  and  T.  G.  Tensfeldt 
Hall-Atwater  Laboratories  of  Chemistry 
Wesleyan  University 
Middletown,  CT  06457 

J.  A.  Montgomery,  Jr. 

United  Technologies  Research  Center 
East  Hartford,  Connecticut  06108 

ABSTRACT 

The  potential  energy  barrier  for  the  isomerization  of  vinylidene,  :C  =  CH2,  to  acetylene, 
HC  s  CH,  has  been  calculated  by  the  quadratic  configuration  interaction  method  using  a 
6s6p3d2f,4s2pld  basis  of  atomic  pair  natural  orbitals  with  extrapolation  to  the  complete  basis  set  limit 

(CBS-QCI/[...3d2f]  APNO  model).  The  calculated  barrier  (AE|  =  2.2  ±0.5  kcal/mol)  and  energy 

change  from  vinylidene  to  acetylene  (AEq  =  -43.9  ±0.5  kcal/mol)  are  in  occeUent  agreement  with 

+ 

recent  experimental  values  (AEe  =  2  kcal/mol,  and  AEo  =  -43  ±2  kcal/mol).  In  spite  of  the  smaU 

barrier  height  for  this  strongly  exothermic  reaction,  the  transition  state  is  located  halfway  between  the 
reactant  and  product  in  apparent  violation  of  the  Hammond  postulate.  The  potential  energy  surface 
for  the  isomerization  can  best  be  understood  by  considering  two  distinct  processes.  The  location  of  the 
transition  state  is  determined  by  the  hydrogen  migration  which  has  a  significant  barrier,  whereas  the 
exothermicity  results  from  the  conversion  of  the  lone  pair  of  vinylidene  to  a  ;r  bond  in  acetylene.  Each 
of  these  processes  individually  satisfies  the  Hammond  postulate. 


Supported  in  part  by  AFOSR  under  Contract  F4962fi-89-C-0019. 
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